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ABSTRACT 



We present a comprehensive study of MHD waves and instabilities in a weakly 
ionised system, such as an interstellar molecular cloud. We determine all the criti- 
cal wavelengths of perturbations across which the sustainable wave modes can change 
2 ' radically (and so can their decay rates), and various instabilities are present or absent. 

Hence, these critical wavelengths are essential for understanding the effects of MHD 

■ waves (or turbulence) on the structure and evolution of molecular clouds. Depending 

on the angle of propagation relative to the zeroth-order magnetic field and the physical 
parameters of a model cloud, there are wavelength ranges in which no wave can be 

^ \ sustained as such. Yet, for other directions of propagation or different properties of a 

■ model cloud, there may always exist some wave mode(s) at all wavelengths (smaller 
CO \ than the size of the model cloud). For a typical model cloud, magnetically -driven 

ambipolar dijfusion leads to removal of any support against gravity that most short- 
wavelength waves (or turbulence) may have had, and gravitationally- driven ambipolar 
CQ ■ diffusion sets in and leads to cloud fragmentation into stellar-size masses, as first sug- 

gested by Mouschovias more than three decades ago - a single-stage fragmentation 
theory of star formation, distinct from the then prevailing hierarchical fragmenta- 
tion picture. The phase velocities, decay times, and eigenvectors (e.g., the densities 
and velocities of neutral particles and the plasma, and the three components of the 
I magnetic field) are determined as functions of the wavelength of the disturbances in 

^ ■ a mathematically transparent way and are explained physically. Comparison of the 

results with those of nonlinear analytical or numerical calculations is also presented 
where appropriate, excellent agreement is found, and confidence in the analytical, 
linear approach is gained to explore phenomena difficult to study through numerical 
simulations. Mode splitting (or bifurcation) and mode merging, which are impossible 
in single-fluid systems for linear perturbations (hence, the term "normal mode" and 
the principle of superposition), occur naturally in multifluid systems (as do transi- 
tions between wave modes without bifurcation) and have profound consequences in 
the evolution of such systems. 
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1 INTRODUCTION - BACKGROUND 

A typical molecular cloud which has not yet given birth to stars is a cold (T ~ 10 K) but complex, partially ionised system, in 
which self-gravitational and magnetic forces are of comparable magnitude, with thermal-pressure forces becoming important at 
high densities ( > 3 x 10* cm~^) or along magnetic field lines. Mouschovias (1976) showed that, barring external disturbances, 
if the magnetic field were to be frozen in the matter, interstellar clouds that have not yet given birth to stars would remain 
in magnetohydrostatic (MHS) equilibrium states. However, ambipolar diffusion (the relative motion of neutral particles and 
charged particles attached to magnetic field lines) is an unavoidable process in partially ionised media. It reveals itself in 
two distinct ways, depending on whether it is magnetically or gravitationally driven (see discussion in § 4). The two kinds of 
ambipolar diffusion acting together initiate fragmentation and star formation in molecular clouds (Mouschovias 1987a). 

In this fragmentation theory, the evolutionary (or fragmentation, or core formation) timescale is the gravitationally- 
driven ambipolar-diffusion timescale, t^d. This does not mean, however, that it takes a time equal to r^^ to form stars. 
The star- formation timescale can be a fraction or a multiple of r^^ , depending on the mass-to-flux ratio of the parent cloud 
and the degree to which hydromagnetic waves (HM) contribute to the support of the cloud: the closer to its critical value 
the mass-to-fiux ratio is and/or the greater the contribution of HM waves to cloud support, the faster the evolution and the 
shorter the star-formation timescale (e.g., see Mouschovias 1987a; Fiedler & Mouschovias 1993, Fig. 9a; Ciolek & Basu 2001; 
Tassis & Mouschovias 2004, Fig. 4).|j The collapse retardation factor, v^^ = Tjf/rni (where % is the free-fall time and Tni the 
neutral-ion collision time), is the factor by which magnetic forces slow down the contraction relative to free-fall. This was a new 
theory of fragmentation (or core formation), initiated by the decay, due to ambipolar diffusion, of relatively small- wavelength 
perturbations (Mouschovias 1987a, 1991a). 

Magnetic braking operates on a timescale shorter than the ambipolar-diffusion timescale and even the free-fall timescale, 
and keeps a cloud (or fragment) essentially corotating with the background up to densities ~ 10* — 10® cm~"^ and thus resolves 
the angular momentum problem of star formation. More specifically, the entire range of periods of binary stars from 10 hr to 
100 yr was shown to be accounted for by this self-initiated mode of star formation (Mouschovias 1977). Even single stars and 
planetary systems become dynamically possible (Mouschovias 1978, 1983). 

Star formation, whether self- initiated or triggerred (e.g., by a spiral density wave, or the expansion of an Hll region 
or a supernova remnant; see review by Woodward 1978) is an inherently nonlinear process. Ambipolar-diffusion— initiated 
star formation has been studied analytically (Mouschovias 1979, 1991a, b) and numerically using adaptive grid techniques in 
axisymmetric geometry up to densities ~ 10^*^ cm~^, by which isothermality begins to break down (Fiedler & Mouschovias 
1992, 1993; Ciolek & Mouschovias 1993, 1994, 1995; Basu & Mouschovias 1994, 1995a, b). More recently, these calculations 
were extended into the opaque phases of star formation (Desch & Mouschovias 2001; Tassis & Mouschovias 2007a, b, c; Kunz 
& Mouschovias 2009, 2010). The key conclusions of the earlier analytical calculations have been verified and numerous new, 
specific, quantitative predictions have been made, many of which have been confirmed by observations (e.g., see Crutcher et 
al. 1994; Ciolek & Basu 2000; Chiang et al. 2008; reviews by Mouschovias 1995, 1996). One may nevertheless take a step 
back from the relatively complicated numerical calculations and ask which results, if any, of the nonlinear simulations can 
be recovered with a linear analysis. With one's confidence increased in the validity of the linear approach (within some self- 
evident limits), one may then make new predictions concerning phenomena that cannot be or have not been included yet in 
the nonlinear calculations. 

The propagation, dissipation, and growth of perturbations in a physical system depends both on the nature of the 
perturbations and the properties of the system. Hydromagnetic waves (or MHD turbulence) seem to play a significant role 
in molecular clouds on lengthscales typically greater than ~ 0.1 pc. They have been shown to account quantitatively for the 
observed supersonic but subAlfvenic spectral linewidths: an observational almost-scatter diagram of linewidth versus size is 
converted into an almost perfect straight line if plotted in accordance with a theoretical prediction by Mouschovias (1987a), 
which relates the linewidth, the size, and the magnetic field strength of an observed object (see Mouschovias & Psaltis 1995, 
Figs. 1 and 2, and update by Mouschovias et al. 2006, Figs. 1 and 2). 

Using a linear analysis as a first step in understanding nonlinear phenomena is not, of course, a new idea. Jeans (1928) 

^ A number of authors assume that molecular clouds are highly magnetically supercritical (e.g., Mac Low & Klessen 2004 and references 
therein; Lunttila et al. 2008, 2009). Zeeman observations (e.g., Crutcher 1999) are sometimes used as the observational justification of 
that assumption. However, geometrical corrections are ignored in those assumptions. The corrections are necessary because (1) only the 
line-of-sight component of the magnetic field is measured, and (2) the measured column density of a cloud flattened along the magnetic 
field lines is statistically greater than that needed for the calculation of the cloud's mass-to-flux ratio (e.g., see Shu et al. 1999). However, 
even if the Zeeman observations were to be taken at face value, without any geometrical correction at all, they still do not reveal highly 
supercritical molecular clouds. In order to obtain a magnetically supercritical molecular cloud model of mass ~ 1840 Mq and mean 
density ~ 100 cm~^, Lunttila et al. assume a magnetic field of only 0.69 /iG or 0.34 /^G. These values are smaller by a factor of ~ 10 
than even the observed strength of the magnetic field in the general interstellar medium 6 /^G (Heiles & Crutcher 2005) . There is no 
conceivable physical mechanism that can possibly increase the density of a forming cloud by 2 — 4 orders of magnitude while decreasing 
its magnetic field strength, especially by that large a factor (~ 10). Hence, the main assumption of many turbulence simulations, i.e., 
that molecular clouds are highly supercritical, has neither an observational basis nor any theoretical justification. 
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used it to obtain his famous instability criterion for the collapse of a cloud against thermal-pressure forces. Hardly any 
astrophysical system exists whose stability with respect to small-amplitude disturbances has not yet been studied by using 
at least an idealized, mathematically tractable model of the physical system. A magnetically supported molecular cloud, 
however, defies a simple linear analysis. First, no realistic equilibrium states have been obtained by analytical means. Second, 
to study the role of ambipolar diffusion in star formation, one must use at least the two-fluid magnetohydrodynamic (MHD) 
equations governing the motions of the neutral particles and the plasma (ions and electrons). In fact, as shown by Ciolek 
& Mouschovias (1993), charged (and even neutral) grains play a very significant role in the ambipolar-diffusion— initiated 
protostar formation. One then has to use at least the four-fluid (neutral molecules, plasma, negatively-charged and neutral 
grains) MHD equations even for a linear analysis to be realistic and relevant to typical molecular clouds. In this paper we use 
the two-fluid MHD equations to study the propagation, dissipation, and growth of HM waves in an idealized model molecular 
cloud. In a subsequent paper we consider the effects of the grain fluid(s). 

Langer (1978) studied the stability of a model molecular cloud (infinite in extent and uniform in density and magnetic 
field) with respect to small-amplitude, adiabatic perturbations in the presence of ambipolar diffusion. For propagation along 
the magnetic field lines, he recovered, as one would expect intuitively, the Jeans dispersion relation and instability in the 
absence of the magnetic field. He then investigated the wave propagation perpendicular to the field lines. He showed that the 
Jeans instability is still present, that the critical wavenumber for instability is independent of the magnetic field strength, 
but that the growth rate depends on both the field strength and the degree of ionisation. Aside from two spurious curves in 
his Figure 1, which exhibits the growth rate and decay time of some modes, our results for propagation of the low- frequency 
modes perpendicular to the magnetic field are in agreement with Langer's - he ignored the high-frequency ion modes. Yet 
even in this, previously studied case, we offer new analytical expressions and new physical insight and interpretation of the 
results. Moreover, we present not only the eigenvalues (frequencies or, equivalently, phase velocities) as functions of wavelength 
but the eigenvectors as well (i.e., material velocities, densities, magnetic-field components, etc.). We also study propagation 
at arbitrary angles with respect to the magnetic field, and we offer a thorough discussion of the wave modes, not just the 
ambipolar-diffusion-induced instability. 

Pudritz (1990) revisited Langer's problem (with the minor difference of considering isothermal perturbations) but intro- 
duced a new effect: he assumed that there exists a power-law spectrum of small-amplitude waves, and then he studied the 
effect that this spectrum has on the ambipolar-diffusion-induced, Jeans-like instability. He concluded that the slope of the 
spectrum (considered as a function of wavelength) has an important effect on the growth rate of the instability; the steeper 
the spectrum, the greater the growth rate. (The growth rate of gravitationally-driven ambipolar diffusion, however, cannot 
possibly exceed the free-fall rate.) 

Several other papers have appeared in print since 1990, studying different aspects of weakly ionised systems, focusing 
usually on the stability of certain MHD modes or shocks, especially as it may relate to the formation of structures in molecular 
clouds and/or on the effect of the grain fluid(s) on the allowable wave modes or shocks (e.g., Wardle 1990; Balsara 1996; Zweibel 
1998; Kamaya & Nishi 1998, 2000; Mamun & Shukla 2001; Cramer et al. 2001; Falle & Hartquist 2002; Tytarenko et al. 2002; 
Zweibel 2002; Ciolek et al. 2004; Lim et al. 2005; Oishi & Mac Low 2006; Roberge & Ciolek 2007; van Loo et al. 2008; Li & 
Houde 2008). In this paper we present a general theory of the propagation, dissipation and growth of MHD waves in partially 
ionised media in three dimensions, with emphasis on mathematical transparency of the formulation and analytical solution 
of the problem, the physical understanding and interpretation of all modes, including their eigenvectors, the many critical 
wavelengths that exist and which separate regimes dominated by different waves or instabilities, and on specific features 
relevant to the evolution of molecular clouds. As mentioned above, even when a particular result agrees with previous work, 
we offer new insight into its physical understanding. 

In § 2 we present the equations governing the behaviour of a weakly ionised, magnetic, self-gravitating interstellar 
cloud. The equations are linearised, Fourier-analyzed, and put in dimensionless form. The free parameters of the problem 
are identified, their physical meaning explained, and their typical values given. The different hydromagnetic modes and their 
dependence on wavelength for different directions of propagation relative to the magnetic field are calculated and explained 
physically in § 3. Analytical expressions for the phase velocities, damping timescales, growth timescales, including critical 
or cutoff wavelengths, are also obtained. A physical discussion of the eigenvectors is an integral part of this presentation. 
Section 4 summarizes some of the results and their relevance to the formation of protostellar fragments (or cores) and to 
other observable phenomena. It also gives in two Tables all the critical wavelengths and the ranges of wavelengths in which 
different modes can exist in molecular clouds, for propagation parallel, perpendicular, and at arbitrary angles with respect to 
the magnetic field. 



The plasma force equation in Pudritz's (1990) paper (eqs. [2.2b], [3.4], [A4], and [AlO]) and in Langer's (1978) paper, eq. (6), contains 
an error. The thermal-pressure force should be multiplied by a factor of 2 to account for the presence of electrons. Although this omission 
does not affect Langer's results because he did not consider the ion modes, it does introduce errors in some of the ion modes considered 
by Pudritz. 
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2 FORMULATION OF THE PROBLEM 
2.1 Basic Equations 

We consider a weakly ionised medium (e.g., an interstellar molecular cloud) consisting of neutral particles (H2 with a 20% 
helium abundance by number; subscript n), electrons, and singly-charged positive ions (subscript i). For specificity we assume 
that the ions are molecular ions (such as HCG"''); for the densities of interest in this paper (~ 10^ cm~^), this is sufficient since 
atomic ions (such as Na+ or Mg+) are less abundant (a result of depletion of metals in dense clouds) and, in any case, they 
have masses comparable to that of HCO^ (for more detailed treatments of the chemistry, see Ciolek & Mouschovias 1995, 
1998, or the appendix of Mouschovias & Ciolek 1999). Interstellar grains, which have been shown to have significant effects 
on the formation and contraction of protostellar cores (Ciolek & Mouschovias 1993, 1994, 1995) and in the opaque phase of 
star formation (Tassis & Mouschovias 2007a, b, c; Kunz & Mouschovias 2009, 2010) are neglected in this analysis; they are 
accounted for in a subsequent paper. 

The magnetohydrodynamic (MHD) equations governing the evolution of the above two-fluid system are 



^+V-(p„«„) 

|i + v.(p..o 



= 0, 



^ Pn Odr 2 
'"iCcn Pi ' 



Pn 



d_' 

dt 

d_ 

dt 



T- T 



Pn 

I 

dB 

'dt 



V B = 



7 

(Vx-B) xB _ pi_ 

AtV Tin 



Vx (viXB) 

kBT 

Pn , 

pmn 
inGpn, 
0, 



pn 



{Vi - Vn) ; 



where pa and Va are, respectively, the density and velocity of species a, 
d' 



dt 



d 



(la) 
(lb) 

(Ic) 

(Id) 

(le) 

(If) 

(Ig) 
(Ih) 
(li) 

(2) 



is the time-derivative comoving with species a, %p the gravitational potential, Pn the neutral pressure, B the magnetic field, 
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the neutral-ion and ion-neutral mean collision (i.e., momentum-exchange) times. The quantity G is the universal gravitational 
constant, and fca is Boltzmann's constant; is the mean mass of a neutral particle in units of the atomic-hydrogen mass 
and is equal to 2.33 for a H2 gas with a 20% helium abundance by number. The quantities C,^-^ and Qdr in the ion mass 
continuity equation \Vo\ are, respectively, the cosmic-ray ionisation rate and the coefficient for dissociative recombination of 



molecular ions and electrons (in cm 



In writing equation (|lb[l . we have used the condition of local charge neutrality 



e(ni — rie) = 0, where rii and rie are, respectively, the number densities of ions and electrons. We may use the assumption of 
local charge neutrality because the various HM modes of interest here have frequencies much smaller than the electron plasma 
frequency ajp,c = (4-7rnce^/mc)^''^ = {A-nx^nne^ /m^Y^'^ = 5.64 x 10^ (a::c/10"'')^''^(nn/10^ cm"^)^''^ s"^ (where = n^/un is 
the abundance of electrons relative to the neutrals, and is equal to the degree of ionisation for weakly-ionised systems); hence, 
any excess charge density is quickly shielded by the mobile electrons, so that no — ni for timescales >ct;~g. Since we neglect 
the effects of grains in this paper, capture of ions onto grains is not included on the right-hand side of equation (|lb|) as a sink 
term for ions. 



Heat conduction and viscosity are not important for the densities and lengthscales of interest 



10^ - 10" 



L ~ 10 ^ — 10^ pc), and are therefore ignored as possible sources of heating/cooling in the model clouds. (Note that the 
left-hand side of equation [le] is equal to pnT[d/ dt]nS , where S is the entropy per gram of matter.) 
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Because pn ^ Pe, Pi, we include only the neutral density as a source term in Poisson's equation (|llip . Similarly, the 
gravitational and thermal-pressure forces (per unit volume) on the plasma (ions and electrons) have been neglected in the 
plasma force equation Hld[l . One can easily show that, for the physical conditions in typical molecular clouds, they are 
completely negligible in comparison to the magnetic force exerted on the plasma, except in a direction almost exactly parallel 
to the magnetic field. Ignoring them parallel to the magnetic field implies that we are neglecting the ion acoustic waves and 
the (extremely long-wavelength) Jeans instability in the ions. 

The quantity {aw}^-^_^ in equations (|3a|l and (|3b|l is the elastic collision rate between ions and neutrals. For HCO^ — H2 
collisions, {aw)^j^^ = 1.69 x lO^'"* cm'' s~^ (McDaniel & Mason 1973). The factor 1.4 in equations (|3a|) and (|3b|l accounts for 
the inertial effect of He on the motion of the neutrals (for a discussion, see § 2.1 of Mouschovias & Ciolek 1999). 



2.2 Linear System 

To investigate the propagation, dissipation, and growth of HM waves in molecular clouds, we follow the original analysis by 
Jeans (1928; see also Spitzer 1978, § 13.3a; and Binney & Tremaine 1987, § 5.1), and assume that the zeroth-order state is 
uniform, static (i.e., Va = 0), and in equilibrium. [f| We consider only adiabatic perturbations; therefore, the net heating rate 
(Fn — An) on the right-hand side of equation (|le|l vanishes. 

We write any scalar quantity or component of a vector gtot(f',i) in the form gtot(T',t) = go + QiXyt), where 50 refers to 
the zeroth-order state, and the first-order quantity q satisfies the condition |g| <C |g()|. We thus obtain from equations (|la|) - 
(fll)) the linearised system 
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Equation H4b|l has been simplified by using the relation 

2 

" /irriH V nii 

which expresses equilibrium of the ion density in the zeroth-order state, as a result of balance between the rate of creation 
of ions from ionisation of neutral matter by high-energy {E > 100 MeV) cosmic rays and the rate of destruction of ions by 
electron— molecular-ion dissociative recombinations. This relation allowed us to replace (^^^ by adiipi,o/mi)^ {fimn/ pn.o) in 
equation (|4b[) . The quantity Xi.o = rufi/rin.o in equation (|4b|) is the degree of ionisation (where ni,o and nn,o are the number 
densities of ions and neutrals in the unperturbed state). For an ideal gas (with only translational degrees of freedom), 7 = 5/3 
in equation (|4e|) . 

We seek plane- wave solutions of the form q{r, t) — g exp(ife ■ r — itot), where k is the propagation vector, u the frequency. 



^ It is well known that the assumption that the gravitational potential is uniform in the zeroth-order state is not consistent with Poisson's 
equation l llhll (e.g., see Spitzer 1978; Binney & Tremaine 1987). However, it is not well known that, for an infinite uniform system, no such 
inconsistency exists; i.e., there is no net gravitational force on any fluid element, hence this state is a true, albeit unstable, equilibrium 
state. 
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Figure 1. Coordinate system used in analysing the linearised hydromagnetic equations. The x-axis is aligned with the propagation 
vector k, and the magnetic field Bq is in the (a;, 2:)-plane at an angle 9 with respect to k (see eqs. |9a| - |9c|). 



and q the amplitude (in general, complex) of the perturbation. Equations (|4a|) - (l4i| reduce to 
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a,o = ^ = (7) 
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is the adiabatic speed of sound in the neutrals, and 

r,,^„ = (4^Gpn,o)-'/' (8) 

is the (one-dimensional) neutral free-fall timescale. The quantities T, Pn, and have been eliminated by using equations l|4ep . 
(Hgl, and (|4hp . respectively. 



2.3 The Dimensionless Problem 

We put equations (|6a[) - (|6f[l in dimensionless form by adopting Pn,o, -Bo, Tjj q, and Ca.o as units of density, magnetic field 
strength, time, and speed, respectively. The implied unit of length is Ca,o''"tt o' which is proportional to the one-dimensional 
thermal Jeans lengthscale Aj^th (see § 3.1.1). For convenience, we adopt a cartesian coordinate system such that the propagation 
vector k is in the a;-direction and the zeroth-order magnetic field Bq is in the {x, 2:)-plane at an angle 6 with respect to k (see 
Fig. 1). Then the three unit vectors are 

fe 

= -, (9a) 
- ^"^^ (9b) 
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One may write Bq — e^Bo cos9 + BzBo sind, and the dimensionless form of equations (|6al) - (|6f|) can be written in component 
form as 

Cjpn = kVn.x, (10a) 
liPi = Pi.ofeSi.a: + i/5f_o3m,drPn — i2pi,0«m,drPi, (10b) 
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Note that equation (jlOll) is redundant in that it gives the same information as equation (|10ep . namely, that there cannot be 
a nonvanishing component of the perturbed magnetic field in the direction of propagation. 
The dimensionless free parameters appearing in equations (|10a|l - (|10ip are given by 

Tni,o = = 0.506 , (11a) 

n., = ^ = 6.29xlO-fi5^V'\ (lib) 

- "dr 1 /ll ^ inS "dr W Wn.O \l/2/29amu\ 
Qm.dr = Pn,OT-ff o = 1-41 X 10 —r 3— ^ —3 ; (Hd 

rui • V 10 " cm'^ s ^/ V lO'^ cm ■^Z \ mi / 

they represent, respectively, the neutral-ion collision time, the ion-neutral collision time, the ion Alfven speed i'A,i,o = 
Bo/ (inpifi)^^'^ , and the electron— molecular-ion dissociative recombination rate per unit ion mass. In evaluating the nu- 
merical constants in equations (jllal) - (|lld|) we have used fi — 2.33 and 7 — 5/3; we have also normalized the ion mass to 
that of HCO^ (= 29 amu). For any given ion mass mi and mean mass per neutral particle fi (in units of mn), the ion mass 
fraction pi,o and the cosmic-ray ionisation rate ^'''^ not free parameters in the problem; the former is determined by the 
ratio fin.o/'fni.o and the latter by the product pf^o'^m.dr, where 3m, dr is the dimensionless dissociative-recombination coefficient 
for molecular ions. 

We note that fni,o = l/i^ufi, where lygfi is the collapse retardation factor, which is a parameter that measures the effec- 
tiveness with which magnetic forces are transmitted to the neutrals via neutral- ion collisions (Mouschovias 1982), and appears 
naturally in the timescale for the formation of protostellar cores by ambipolar diffusion (e.g., see reviews by Mouschovias 
1987a, § 2.2.5; 1987b, § 3.4; 1991b, § 2.3.1; and discussions in Fiedler & Mouschovias 1992, 1993; Ciolek & Mouschovias 
1993, 1994, 1995; Basu & Mouschovias 1994, 1995). It is essentially the factor by which ambipolar diffusion in a magnetically 
supported cloud retards the formation and contraction of a protostellar fragment (or core) relative to free fall up to the stage 
at which the mass-to-flux ratio exceeds the critical value for collapse. It is discussed further in § 3.2.1. 

Equations (jlOap - (|10k|) govern the behaviour of small-amplitude disturbances in a weakly ionised cloud; they (without eq. 
|10ej l form a 10 x 10 homogeneous system. In general, the dispersion relation a;(fc) can be obtained by setting the determinant 
of the coefficients equal to zero. To each root (eigenvalue) of the dispersion relation there corresponds an eigenvector (or 
"mode"), whose components are the dependent variables appearing in equations (|10ap - (|10k|) . (Note that, once the dependent 
variables in eqs. |10a| - |10k| are known, one may use eqs. [4e], [4g], and [4h| to solve for the perturbed quantities T, Pn, and tp, 
respectively.) Since, in general, uj is complex, modes with Im{Li3} < decay and those with Im{d;} > grow exponentially in 



8 Mouschovias, Ciolek, and Morton 



time. In what follows we investigate the propagation, dissipation, and growth of the allowable HM modes in typical interstellar 
molecular clouds. 

3 SOLUTION, PHYSICAL INTERPRETATION, AND APPLICATIONS 

For specificity, we consider a representative molecular cloud of density rin.o — 2x 10^ cm"'', magnetic field strength Bo ~ 30 fiG, 
temperature T = 10 K, dissociative recombination rate cvdr = 10~® cm'' s~^, and cosmic-ray ionisation rate ^^j^ = 5x 10"^^ s~^, 
implying a degree of ionisation a;i,o — 1.58 x 10~^ (and, hence, ion mass fraction pi,o ~ 1.97 x 10~^). The unit of time r,.f ^ 
(see eq. [8]) for this model is equal to 3.92 x 10^ yr, and the unit of speed is Ca,o ~ 0.243 km s~^ (see eq. [?])• Hence, the 
unit of length is Ca,o''"ff g ~ 9.72 x 10~^ pc. The four (dimensionless) free parameters of the problem (see eqs. |lla| - llldp 
are: the ion-neutral collision time fin.o ~ 4.45 x 10~^, the neutral-ion collision time fni.o = 0.226, the Alfven speed in the ions 
VA,i,o = 2.81 X 10'^, and the dissociative-recombination coefficient Om.dr = 1.99 x 10®. (These imply a dimensionless cosmic-ray 
ionisation rate (^^j^ =6.19 x 10~^.) 

In the following subsections we present the solutions for propagation along {9 = 0°), perpendicular {9 = 90°), and at 
intermediate angles (6 = 45°, 10°, and 80°) with respect to the unperturbed magnetic field Bo (see Fig. 1). We use the 
velocity vector v in relation to k as defining the polarisation of each wave mode. Modes that have only Vn,x, Vi,x 7^ are said 
to be longitudinally polarised; it follows from equation (|10ap that these modes are compressible, i.e., pn 7^ 0. Modes that have 
Vii,x, = are said to be transversely polarised; they are incompressible, i.e., pn = 0. Note in what follows that, since the 
thermal pressure in the plasma has been neglected, no ion sound waves are present. 



3.1 Propagation Along Bq {k \\ Bo) 

For 9 = 0°, equations |10aj - [10k| become uncoupled in the three mutually orthogonal directions Sx, e^, and e^. The modes 
polarised in the x-direction are given by (see eqs. [lOaj - |10d| 'l 

ujpn = kvn,x, (12a) 

ti/Si = pifikVi,x + ipfofim.drPn — i2/5i,0(5:m,drPi , (12b) 

CjV^^x = I — i ) Pn — ^J—V^,x + -J—Vi^x, (12c) 
V kj Tnifi Tni,0 

z i, 

I^Vi^x = Vn,x — Vi,x- (12d) 

Tin.O T-in,0 

Equations (jlOf|l - (I10h|l yield for the modes with motions in the j/-direction 
z i 

^Vn,y = —- V^,y + V\^y, (13a) 

OjVi,y = -r—Vn,y — -r—Vi,y — VA,i,okBy, (13b) 
■7"in,0 Tin,0 

CjBy = —kvi^y, (13c) 
while equations (|10i[l - H10k[l for the modes polarised in the z-direction become 

OJVn,z = -~-r—Vn,z + ^T—Vi^z, (14a) 
Tni.O Tni,0 

<iCi,z = -^-—Vn,z — -r—Vi,z — y%iokBz, (14b) 

T"in,0 Tin.O 

ujBz = —kvi^z- (14c) 



3.1.1 Longitudinal Modes for k \\ Bo: Eigenfrequencies and Eigenvectors 

From equations (|12a|) - ()12d|) the dispersion relation for the longitudinal modes is 



{UJ + 2ipi^oQm,dr 



Ul +t { h 

Tin,0 Tnifi 



I 

T'in.O 



= 0. 



(15) 



The four eigenvalues and eigenvectors as functions of dimensionless wavelength A (= 2Ti/k) are obtained from direct numerical 
solution of the eigensystem (|12ap - (|12dl) and displayed in Figures 2 and 3. Figure 2a shows the absolute value of the phase 
velocity (= Ur/k, where Cjr = Re{a;}). The damping (or dissipation) timescale (= uji < 0, where uji = Im{ci)}) is 

exhibited in Figure 2b; growth times t^^ (= l/Coi, uJi > 0) are shown in Figure 2c. The absolute values of the Wn.i-components 
of the eigenvectors for the various modes are shown in Figure 3a, while those of Vi^x, pn, and pi are displayed in Figures 36, 
3c, and 3d, respectively. Note that, in Figures 3a - 3d, all the eigenvectors have been normalised to unity; i.e., they satisfy 
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Figure 2. Eigenvalues of longitudinal modes as functions of wavelength, normalised to Ca.oTff q = 9.72 X 10~^ pc, at an angle of 
propagation 6 = 0° with respect to Bq. Because of degeneracy, there are a total of four different modes at eacli wavelength. Each curve 
is identified by a label (a mnemonic for the mode it represents) as explained in Table 1. (a) Absolute value of the phase velocity v^, 
normalised to Ca,o = 0.243 km s~^. Also overplotted as boxes with interior x's are values obtained from eq. I I19II . (6) Damping timescale 
T^, in units of ^ = 3.92 X 10^ yr. Also shown (open circles) arc values calculated from eq. I I21II . (c) Growth timescale r^^, normalised 
to Tff (J . Also plotted (open circles) are values calculated using eq. Ij2ip . 
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Figure 3. Magnitudes of eigenvectors of longitudinal rnudes for 9 = 0° as functions of wavelength, normalised as in Fig. 2. As in Fig. 
2, there are a total of four modes. All eigenvectors have been normalised to unity. Each curve is identified by a label (a mnemonic for 
the mode it represents) as explained in Table 1. (a) Longitudinal (a;-) component of the neutral velocity, li^nja:!. (b) Longitudinal (x-) 
component of the ion velocity, \vi^^\. (c) Neutral density |pn|- (d) Ion density |pi|. 
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the condition (|pnp + |pip + iwn.ip + = 1- The first mode for this system of equations is a nonpropagating, ion 

collisional- decay mode, i.e., the ions are streaming through a sea of fixed neutrals. Their motion decays because of ion-neutral 
collisions. It is characterized by v„_^x — Q = pn- Solving equations p2a|) - (|12d|l under these conditions (or, equivalently, solving 
eq. [15] in the limit \u}\ ^ l/fni,o), one finds 



(16) 



Hence, = Q and = fin,o ~ 4.45 x 10~^ (see Fig. 26, line labeled "i,coU"). This mode is independent of wavelength (see 
Figs. 26, 3a — 3d, lines labeled "i,coU"). 

The second mode is one in which density enhancements in the ions rapidly decay by dissociative recombinations of 
molecular ions and electrons; because the degree of ionisation is so small (si.o = 1.58 x 10~^), this mode does not involve any 
motion of the neutrals and leaves the neutral density essentially unchanged (see Figs. 3a and 3c, lines labeled "i,rec"). Solving 
equation (|12b[) with pn = and Vi^x = 0, one finds that = Q and 

= (2(5:m,drPi,0)~^ , (17) 

which is equal to 1.28 x 10"* for the model described here. It is again the case that this mode is independent of k (see Figs. 
26 and 3a - 3d, lines labeled "i,rec"). 

The remaining two longitudinal modes are low-frequency modes, with |a;| <C l/fin,o. Because the inertia of the ions along 
the magnetic field is small {pi.o ^ pn,o), the neutrals are able to sweep up the ions, and, as a result, Vi^x — v-n^x (see Figs. 3a 
and 36). For these conditions, equation (|15p yields the thermal Jeans modes 

(18) 



±k 



1- 



(e.g., see Chandrasekhar 1961, Ch. XIII; or Spitzer 1978, 

1/2 



13.3a). Therefore, for fc > 1, 



1 



(19) 



i.e., the two acoustic waves have the same phase velocity, modified by gravity, but propagate in opposite directions (along 
the field lines). In the limit \ <^ 1, = ±1 and = oo; i.e., these modes are undamped sound waves (recall that the 
unit of speed is Ca,o). At longer wavelengths, gravitational forces become increasingly more important and the phase velocity 
becomes less than unity (see Fig. 2a). The waves are gravitationally suppressed (i.e., = 0) at wavelengths greater than the 
thermal Jeans wavelength 



= 271 



(20) 



(= 27rCa,oT"ff = 0.611 pc, dimensionally) . For A > Aj,th (i.e, k < 1), it follows from equation (11811 that each of the Jeans 
modes splits into two separate, conjugate modes. One is a gravitational growth (or fragmentation) mode, with timescale 

-1/2 

(21) 



1 - 



At 



A 



(see Fig. 2c). This is the classical Jeans instability. As A — >■ cxd, fj,_. 
. The corresponding eigenvector is labeled as "ff+" in Figures 3a 



► 1; dimensionally, this is just the free-fall timescale, 
3d. The other mode is one of exponential decay, with 
damping timescale also given by equation (|2ip (see eq. [18]); it is the curve labeled by "ff— " in Figure 26. The eigenvector, 
also labeled by "ff— ", is shown in Figures 3a - 3d. This mode is one in which an initial density enhancement causes expansive 
motion, opposed by gravity, at such a rate that the density enhancement decreases to zero at the same time that the velocity 
vanishes. Hence, this is a monotonically decaying mode; no wave motion is involved. It is similar to the well-known classical 
cosmological problem of an expanding "flat" universe. Note that, as A — >■ oo, — >• 1. 

In order to better understand the neutral thermal (Jeans) modes, we examine more closely the eigenvectors (and their 
features shown in Figures 3a - 3d). We substitute equation (|18|) in equation (|12a[l . and we use the normalisation condition 
PnP + |'5n,a;P + I'Ci.a;!^ = 1 aud the fact that \vn,x\ = to fiud that 



and 



_P 

3fe2 - 



3fc2 - 2 

From these equations, it follows that 



(22a) 



(22b) 
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Table 1. GLOSSARY OF Labels in Figures 



Label 


Meaning 


acoustic 


neutral acoustic wave 


ijCoU 


ion coUisional-decay mode 


i,rec 


ion dissociative-recombination mode 


fr+ 


Jeans free-fall mode 


ff- 


conjugate Jeans ( "cosmological" ) mode 


i,A 


ion Alfven wave 


n,A 


neutral Alfven wave 


AD 


magnetically-driven ambipolar-diflusion mode 


n,coll 


neutral coUisional-decay mode 


i,ms 


ion magnetosonic wave 


n,ms 


neutral magnetosonic wave 


PD 


neutral pressure-driven diffusion mode 


AD,fr 


neutral gravitationally-driven AD fragmentation mode 


ijfast 


ion fast wave 


n,slow 


neutral slow wave 


n,fast 


neutral fast wave 



(a) as fc — >■ oo (A — >■ 0), |pn| = |^n,a;| ~ \vi.x\ ^ / V^) = 0.577, as seen in Figures 3a - 3c (curves labeled "acoustic"); 

(b) as — >■ 1 (A — >■ 2tt), \pn\ — >■ 1 but \vn,x\ ~ Ivi.xl — >■ 0, as also seen in Figures 3a - 3c; 

(c) for k < (2/3)^''^ (i.e., A > (3/2)^''^ 2n = 7.70), p„ becomes imaginary (its absolute value is shown in Fig. 3c); 

(d) for (2/3)^/^ < fc < 1 (i.e., 27r < A < 7.70), Vn,x = Si.x is imaginary (Figs. 3a and 36 show the absolute values of these 
velocities). 

For the convenience of the reader, Table 1 contains a list of all abbreviations (and their meaning) used to label the curves 
in all the figures of this paper. 



3.1.2 Transverse Modes for k \\ Bq: Eigenfrequencies and Eigenvectors 



Comparing the systems of equations for the modes with motions only in the y and only in the z directions, p3a|) - p3c[) and 
H14ap - (|14c[) . we note that they are identical. Hence, for propagation along the field {6 — 0°), uj{k) is degenerate for these, 
transverse modes. Figure 4a displays the absolute value of their phase velocity, obtained by solving the dispersion relation 

, 2 



~3 I ■ 



1 ^ 1 



VAA,ok 



0. 



(23) 



(Note that, because of the degeneracy, there are four different modes.) Damping timescales are shown in Figure 46 as functions 



of A. None of the modes are unstable. The absolute values of the eigenvectors are exhibited in Figure 5: 



and \vn,z\ in 



Figure 5a, 



and \vi.z\ in Figure 56, and \By\ and \Bz\ in Figures 5c and 5d. Note that the \By\ (and \Bz\) axis in Figure 



5d is logarithmic in order to show the behaviour of the modes at small wavelengths. 

From the dispersion relation (|23|) and Figures 4a and 5a, it is evident that small-wavelength, high-frequency (|a;| > l/fin,o) 
ion modes propagate with Vn,y, Vn,z — 0. Solving equations (I13a|l - p3c|l (or, equivalently, eqs. |14aj - |14c) ) in these limits 
yields 

1/2 

(24) 



uj = ±{;a i ok 



1 



2SA,i,o'nii,ofc 



2fin,0 



(In deriving eq. [24] we have used the fact that l/fin,o ^ l/'fni.o; see eqs. [lla| and [TTb].) For A less than the ion Alfven 
cutoff wavelength 



AaJ = 47r{;A,i,0'?in,0, 

waves propagate with 



V4, = ±WA,i,0 



A 

aTI 



1/2 



(25) 



(26) 



Aa.i ~ 1.57 X 10 ^ (i.e., Aa.i ~ 1.53 x 10 pc) for the model cloud parameters specified at the beginning of § 3. Hence, 
for A <c; Aa,!, the waves are Alfven waves, with — ±WA,i,o = 2.81 x 10"^ (see Fig. 4a, curve labeled "i,A"). There are four 
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Figure 4. Eigenvalues of transverse modes as functions of wavelength, at an angle of propagation 9 = 0° with respect to Bo- All 
normalisations are as in Fig. 2. At each wavelength there are four different modes in all: two with motions in the y direction, and two 
with motions in the z direction, (a) Absolute value of phase velocity, \v^\. Phase velocities resulting from eqs. I I26I I and II33I I are also 
displayed (circles with interior crosses, and circles with interior asterisks, respectively), (b) Damping timescale Tj. Also shown are values 
(downward- facing triangles) calculated from eq. Il29al l. and values (crosses) calculated from eq. 1134a l l. 
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Figure 5. Magnitudes of eigenvectors of transverse modes for 9 = 0° as functions of wavelength. All labels and normalisations are as in 
Fig. 3 and Table 1. As in Fig. 4, there are four modes at each wavelength, (a) Transverse {y- and z-) components of the neutral velocities, 
|Dn,s/|) and It'll, z|- {b) Transverse (j;- and z-) components of the ion velocities, \vi^y\ and \vi^z\- {c, d) Transverse (y- and z-) components 
of the magnetic field, and \Bz\ shown on two different vertical scales, which bring out different important features (see text). 



MHD Waves in Weakly Ionised Media 15 



waves in all. The two polarised in the j/-direction are normal, shear Alfven waves, and the two polarised in the z-direction are 
modified Alfven waves. (At 6* > the latter set of waves are fast waves.) All the waves are damped on the timescale — 2fin,o 
(see eq. [24| and Fig. 4&, curve labeled "i,A") because of collisions with the neutrals. It is noteworthy that the damping time 
is longer by a factor of 2 than that (fin,o) referring to the dissipation (momentum exchange) of ion streaming motion relative 
to the neutrals. This is so because, although a typical ion indeed loses memory of the collective (wave) motion on a timescale 
fin.o, half of the wave energy is stored as potential energy in the magnetic field. Therefore, it takes twice as long for collisions 
to damp the wave than it takes them to damp ion streaming. 

At \ = Aa,! the ion Alfven waves are critically damped. For A > Aa.i, the ion-neutral collision frequency l/h-afi is greater 
than the wave (angular) frequency uj, and the waves can no longer propagate (v^ = 0). At A = Aa.i there is a bifurcation in the 
ion modes (see Fig. 46). Two of the modes (one polarised in the j/-direction and the other in the z-direction), corresponding 
to the negative root in equation (|24p . become ion collisional-decay modes (discussed earlier in § 3.1.1; curves labeled "i,coU" 
in Figs. 4fe and 5a - 5ci), with 

2fin,() 



1- AA,i/A 



1/2 ' 



A > AA,i. 



(27) 



In the limit A 3> Aa.i, By, Bz (see Figs. 5c and 5d), and — )■ Tin,o, just as in equation (|16p . Thus, as A increases, 
magnetic restoring forces on the ions become negligible, and the motion of the ions simply decays on a timescale fin,o because 
of collisions with the neutrals (see Fig. 4b, curve labeled "i,coU"). The remaining two ion modes (again, one polarised in 
the y-direction and the other in the jz-direction) , corresponding to the positive roots of the dispersion relation (eq. [24]), are 
magnetically-driven ambipolar- diffusion modes (see Fig. 46, curve labeled "AD"), in which the ions (and electrons) diffuse 
quasistatically (i.e., with negligible acceleration; this is equivalent to having <^ l/rin.o in eqs. |13b] and [l4b]) relative to 
the stationary neutrals. The damping timescale for these modes is 

2fin,0 



1 - 



1- AA,i/A 



1 1/2 ' 



(28) 



which is the curve labeled as "AD" in Figure 46. In the limit A ^ Aa.i, becomes equal to the ambipolar-diffusion timescale, 
A^ 



47r2I?a,i ■ 



where 



167r2fin.O 



^'A,i,OT"m,0 



(29a) 



(29b) 



is the ion ambipolar-diffusion coefficient. For the values of the free parameters cited at the beginning of § 3, ©a.i = 3.52. 

In Figures 46 and 5 it is also evident that there exists two small-wavelength, low-frequency {\uj\ <C l/fin,o) neutral modes. 
In these modes the plasma and magnetic field lines are essentially stationary (i.e., ~ 0, iii,z ~ and By ~ 0, B^ — 0). 
Under these constraints, equations (|13ap - (|13c|l (and, similarly, eqs. [14aj - |14cj ) yield 



(30) 



These modes are the neutral collisional-decay modes (see curves labeled by "n,coir' in Figs. 4 and 5), in which the motion of 
the neutrals in the y and z directions decays due to collisions with ions that are held fixed in space by the magnetic field; 
= 0, and = fni,o = 0.226 for these modes (see Fig. 46). 
At longer wavelengths (~ Aa.h, see eq. [32] below), collisions between the neutrals and the ions cause the two fluids to 
begin to move together. Hence, magnetic forces on the ions are more readily transmitted to the neutrals. For large enough A, 
the neutrals can sustain a HM wave. At the point A = Aa.ii, the ion ambipolar-diffusion modes and the neutral collisional-decay 
modes merge (see Fig. 46); waves can propagate at longer wavelengths. In the limit l/fin,o, the dispersion relation (|23[l 

has the solution 

1/2 



i^A.i.O 



±VA,n,ok 



Tin,0 
fni,0 



1/2 



k 



^«A,i,OTm,ofc) 



4(rin,o/rni,o) 



* -2 - 72 
2«A,i,0T"in,0K , 



1 — J (^■CA,n,OT'ni,ofc 



1/2 



J -2 - ra 

2"A,n,0'rni,0fe , 



(31a) 



(31b) 



where 'yA,n,o [= '^A,n,o/C'a,o, WA.n.o = Bo/ (47rpn,o)^''^] is the dimensionless Alfven speed in the neutrals. In equation (|31b|) we 
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have eliminated fin,o and VA,i.o in favor of Tni^o and ■yA,n,o (see eqs. [3a] and [3b| . and recall that pn.o = /imHKn.o)- It follows 
from equation H31b|l that > for all A > AA,n, where 

/f- \^^^ 

Aa.ii = 71" I i5A,i,0^ni,0 = •iT"i}A,n,0'?ni,0 (32) 

is the neutral Alfven cutoff wavelength. (This is referred to simply as the Alfven lengthscale in Mouschovias 1987a, 1991a.) 
In the typical model cloud, Aa.h = 2.80, which means that the dimensional Alfven cutoff wavelength is Aa.h = 0.273 pc. For 
A > Aa.ii, equations l|31b|) and (f32)) yield 

2-1 1/2 

(33) 



Vtf, = ±VA,n,0 



Aa,i 



In the limit A ^ Aa.ii, V4, — >■ ±^A,n,o (= ±3.94, since |t;A,n,o| = 0.957 km s~^); these modes are Alfven waves in the neutrals 
(see curves labeled "n,A" in Figs. 4 and 5). The two (oppositely propagating) waves polarised in the y-direction are normal 
Alfven waves. The two waves polarised in the ^-direction are modified Alfven waves; at S > they are fast waves. As seen 
from equation (|31b|) . these modes damp on the ambipolar-diffusion timescale 



= 2r, = 2-^^, (34a) 
where 

A^ 

^'a.n = = i'A.n.oT'iii.O (34b) 



is the neutral ambipolar-diffusion coefficient. Comparing equations (|29b|) and (I34bl) . we note that Da.n ~ {— 3.52 for 
this typical model cloud). We may therefore denote the ambipolar-diffusion coefficient simply as 2?a, without the subscript 
i or n. One should bear in mind, however, that the expression for for the neutral Alfven waves contains an extra factor 
of 2 compared to the expression for for the ion ambipolar-diffusion mode (compare eqs. [34aj and [29a], and make use 
of eqs. [34b] and [29b]). As explained in the case of the ion Alfven waves, it takes twice as long to damp a wave than it 
takes to damp streaming motion (or diffusion). The existence of Aa.h = THA.nTni for Alfven waves in the neutrals was first 
shown by Kulsrud & Pearce (1969), who studied the excitation and propagation of HM waves in the intercloud medium due 
to cosmic-ray streaming (see, also, Parker 1967). Mouschovias (1987a, 1991a, b) discussed the importance of the lengthscale 
Aa.ii and the thermal Jeans lengthscale Aj,th in the formation of protostellar fragments (or cores) in self-gravitating molecular 
clouds. He proposed that fragmentation is initiated by the decay of HM waves due to magnetically-driven ambipolar diffusion 
and the almost simultaneous onset of a Jeans-like instability, due to gravitationally-driven ambipolar diffusion (see discussion 
in § 4 below). 



3.1.3 Transverse Modes for k \\ Bq: Further Discussion of Eigenvectors 

More insight in the physics of the transverse modes can be gained by understanding analytically certain key features of the 
eigenvectors shown in Figures 5a - 5d. In the case of the ion Alfven mode, we may ignore the motion of and the dissipation 
due to the neutrals at short wavelengths and we may use equations (|13b|) . (|13cp . and the dispersion relation uj ~ ^VA,i,ok (see 
eq. [23]), to find that 

^ ^ T^A.i.o. (35) 

By 

(The =p sign on the right-hand side of eq. [35] refers to propagation in the zbs-direction.) Since these modes are incompressible 
(i.e., pn = = Pi), we may use the normalisation condition 

\vi,yf + \Byf = 1 (36) 
to find that 

VA,l.O . V 

Vi.y - T7^ - 1, (37a) 

(«i.i,o + i)' 

and 

By~T -T-r^. (37b) 

('^i.i.o + l)''' 

Since ua.i.o = 2.81 x 10^, it follows from the last equation that \By\ ~ 3.56 x lO""*. This is in agreement with the curve 
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labeled "i,A" in Figure 5d, which is the same as Figure 5c but with a logarithmic scale for \By\ so as to show this small but 
finite value of \By\ at small A. Similarly, the result Vi^y ~ 1 is as shown in Figure 5b, curve labeled "i,A". 

As A increases, remains large and \Vn,y\ (and \By\) small even for A > Aa.i (see Figs. 5a - 5d, curves labeled "i,A"); 
the ion ambipolar-diffusion mode maintains a significant \vi,y\ although ion Alfven waves do not exist for A > Aa.i- As A 
increases toward the neutral Alfven cutoff wavelength Aa.h (see eq. [32]), \vi,y\ decreases and \vn,y\ increases because the ions 
begin to couple to (and induce motions in) the neutrals via collisions. At exactly A — Xa.u, the two velocities become equal. 
For A > Aa.ii, Alfven waves can be sustained by the neutrals, and the magnitudes of all three quantities Vi,y, Vn.y, and By 
are significant. In the long-wavelength limit, the ions are well coupled to the neutrals {vi,y — Vn,y)- Using the normalisation 
condition 

\ijn,yf + \K,yf + \Byf = 2\Vn,yf + \Byf = 1 (38) 



and equations ()13b|) - (|13c|) . we now find that 



= T^A,n,o, (39a) 

By 

^ ""^'y ^ .1/2 (39b) 



and 



(1 + 2Si,„,o; 

By = T- 1 ,1,, . (39c) 

Since for our typical model cloud 'L'A,n,o ~ 3.94, it follows that Vn,y ~ Vi,y ~ 0.696 and By = 0.177, in agreement with Figures 
5a - 5d, curves labeled "n, A" . 

As A increases across Aa.i, the ion Alfven mode disappears (damps) and bifurcates into the (ion) ambipolar-diffusion 
mode and the ion coUisional-decay mode, as discussed in relation to Figure 46; hence the placement of the labels "i,A", "AD", 
and "i,coU" on these curves in Figures 5a - 5d. 

It is also clear from Figures 5a - 5d that the eigenvectors for both the neutral and the ion coUisional-decay modes behave 
exactly as expected on the basis of our discussion of these modes in relation to Figure 46. 



3.2 Propagation Perpendicular to Bo (k _L Bo) 



For 9 = 90°, equations (|10a|l - ([lOkp again decouple into three independent subsystems. The subsystem involving material 
motions in the a;-direction is 

Ljpn = kVn,x, (4Qa) 

<i/5i = Pi,ofc«i,i. + ip?,o"m,drPn - j2pi,0am,drPi, (40b) 

~ l\ i i . 

k - - ] pn~ Vn,x + Vi^x, (40c) 

kj Tnifi Tni.O 

C0Vi,x = -r^Vn.x — -T^Vi.x + v\,iokB;,, (40d) 

Tin,0 ■ ''"in,0 

ujBz = kvs_,x- (40e) 
Motions in the j/-direction are governed by 

ioVa,y = — i5n,H + -r— Si,!/, (41a) 

Tni.O Tni.O 
% Z 

ijJVi,y ~ Vn,y — Vi,y, (41b) 

clBy = 0. (41c) 

Similarly, the equations for the a-components of the neutral and ion velocities are 

i z 

0jVn,z = Vn,z + Wi.z, (42a) 

Tni.O Tni.O 

i z 

OjVi,z = Wn,2 — Wi.z. (42b) 

Tin,0 Tin.O 
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3.2.1 Longitudinal Modes for k _L Bq: Eigenfrequencies 

The dispersion relation is easily obtained from equations (|40a|) 

{UJ + 2i/5i,o(5fm,dr) < 



V-Tin.O Tni.O / I i 



1 -2 72 , 1 / I z 

I'A.i.oK + I k 

Tni.O Tin.O 



d;+ (fc^ - 1) =0. 



(43) 



Because it is a fifth-order polynomial, there are five longitudinal modes in all. The phase velocities, damping timescales, and 
growth timescales for the longitudinal modes are displayed in Figures 6a, 66, and 6c, respectively; eigenvectors are shown as 
functions of A in Figure 7a - 7e. 

At small wavelengths there are again two high-frequency (|a3| ~ l/Tin,o) ion modes, with Vn,x — (they are degenerate 
with respect to the direction of propagation). In this limit, the solution of the dispersion relation is identical with that given 
by equation (|24l) . The modes it represents in this case are ion magnetosonic modes. Because the speed of sound of the ions 
(< Ca.o, since rm > fimu) is negligible compared to the ion Alfven speed (see eq. [llcj l. the phase velocity of the waves for 
A < Aa.i is again given by equation (|26|) and the characteristic decay time is — 2fin,o (see Figs. 6a and 66, curves labeled 
"i,ms"). For A > Aa.i (= 1.57 x 10^^ for the typical model cloud) the ion magnetosonic waves cannot propagate because of 
frequent ion- neutral collisions. Instead, each mode bifurcates (see Fig. 66), just as in the case of ion Alfven waves described in 
§ 3.1.2. One of the two resulting modes is an ion coUisional decay mode (f^ = fin.o; curve labeled "i,coU" in Fig. 66) and the 
other is a magnetically-driven ion ambipolar-diffusion mode, curve labeled "AD" in Figure 66 (f^ = X^/4n'^t>^ for A > Aa.i; 
see discussion preceding eq. [29b] '). 

The third mode is an ion dissociative-recombination decay mode, as discussed in § 3.1.1, in which density enhancements 
in the ions decay rapidly [fj = (2pi,o3m,dr)^^ = 1-28 x 10"'' for the typical model cloud parameters] because of dissociative 
recombinations between molecular ions and electrons (see Fig. 66, curve labeled "i,rec"). The only nonvanishing component 
of the eigenvector for this mode is pi,o (see Figs. 7e and 7a - 7d). 

There also exist two small-wavelength, low-frequency (jci;| <^ l/nn.o) wave modes in the neutrals. For these modes, Vi,x — 
and Bz — (curves labeled "acoustic" in Figs. 7b and 7c), i.e., the neutrals oscillate in the a;-direction in an effectively fixed 
background of ions and magnetic field. Under these circumstances, the solution of the dispersion relation (|43[) is 



Cj — zbfc 









1 - 













1/2 



2rni,o 



The phase velocity can be written as 

1/2 

/AN 

1 



= ± 

where 
A.. 



A 

As. 11 



47rfni,o 



for A < As 



[1 + (2fni,o)2]^/2 



(44) 



(45) 



(46) 



is the acoustic-wave cutoff wavelength; As,n = 2.59 in the typical model cloud. In dimensional form, As.n — 47rCa,oT"ni,o/[l + 
(2fni,o)^]^'''^ ~ 0.252 pc. (Note that for 2rni,o ^ 1, As,n — 47rfni,o; in this limit As.n — 47rCa,oTiii,o = 0.277 pc.) These are sound 
waves in the neutrals modified by gravity and neutral-ion collisions. For A <C As,n, v<t, ~ ±1 (= ±Ca,o, dimensionally) , hence, 
the waves are pure sound waves (see Fig. 6a, curve labeled "acoustic"). They decay on a timescale = 2fiii,o (= 0.452; see 
curve labeled "acoustic" in Fig. 66). For A > As,n, the waves are damped because of collisions with the ions, and = 0, 
much like the damping of the ion Alfven and magnetosonic modes for A > Aa.i. For A > As,n, this modified neutral sound 
wave bifurcates (see Fig. 66) into a pressure-driven diffusion mode (curve labeled "PD") and a coUisional- decay mode (curve 
labeled "n,coir'). We examine these modes in that order. 

The pressure- driven diffusion mode corresponds to the positive root of equation (I44p . The neutrals are diffusing quasistat- 
ically (i.e., with ^ l/rni.o) through a background of effectively stationary ions and magnetic field. The eigenvector for this 
mode is labeled "PD" in Figs. 7o - 7e. For As,n < A < Aj,th, the damping timescale is 



2fni.l 



1 - 4f 2 



(Aj,th/A)' 



- 1 



1/2 ' 



(A/A.J 



Op 



1- A/Aj,th 



(47a) 
(47b) 
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a,0 ff, 



Figure 6. Eigenvalues of longitudinal modes as functions of wavelength, at an angle of propagation 9 = 90° with respect to Bq- All 
normalisations are as in Fig. 2. There are five difl'erent modes at each wavelength, (a) Absolute value of phase velocity, \v^\. Phase speeds 
derived from eq. II26I I are displayed as circles with interior crosses, while speeds obtained from eq. 1145 H are displayed as diamonds with 
interior black circles. Phase speeds determined from eq. II55I I are displayed as asterisks. (6) Damping timescales t^. Also shown are values 
(downward- facing triangles) calculated from eq. I l29al l. as well as values (plus signs) from eq. I l34al l. Results (five-pointed stars) from eq. 
H47ap are also depicted, and so are those (boxes) calculated by using eq. JSGj. (c) Growth timescale t^^. Values predicted (six-pointed 
stars) from the ambipolar-diffusion fragmentation timescale (eq. 1491 ) are also displayed, along with values (boxes) calculated from the 
magnetic Jeans instability mode, eq. H57|l . 
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Figure 7. Magnitudes of eigenvectors of longitudinal modes for = 90° as functions of wavelength. All labels and normalisations are as 
in Fig. 3 and Tabic 1. As in Fig. 6, there are five modes, (a) Longitudinal (a;-) component of the neutral velocity, liinjo;!- {b) Longitudinal 
(x-) component of the ion velocity, \vi^x\- (c) Magnetic field component \Bz\. (d) Neutral density |pn|. 
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Figure 7. Cont. (e) Ion density \p\\. 

where 

= fni,0 (48) 

(= Ca,o''"ni,o dimensionally) is the neutral pressure-driven diffusion coefficient. In obtaining equation (|47b|l we have used 
the fact that f g ^ 1. For the representative model cloud used in this paper, Vp — 0.226 (~ 2?a/16). We note that, at 
A = Aj.th, f J = oo (see Fig. 66, curve labeled "PD"). At this wavelength, the restoring pressure forces in this mode are 
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exactly balanced by self-gravitational forces, and the system is on the verge of gravitational instability {ui = 0). Thus the 
Jeans instability manifests itself at Aj_th even in the presence of a magnetic field, as originally recognized by Langer (1979). 
Ambipolar diffusion allows this to happen, but, because of neutral-ion collisions, the growth time of the instability is longer 
than that of the nonmagnetic Jeans instability (compare eq. [49] below with eq. [21]). For A > Aj,th, density perturbations in 
the neutrals grow exponentially in time as a result of gravitational contraction of neutrals through essentially stationary ions 
attached to rigid magnetic field lines. The growth time for this instability can be obtained easily from equation (|44p under 



the conditions A > Ay th and 4f„ 



1 



Aj,th/A 



< 1: 



(l/rni,0) 



t'ff.O 



(Aj,th/A' ' 



-N 2 ' 

A,j,th/A) 



(49) 



where z^h.o — l/rni,o {— % p/rni.o = 1/0.226 — 4.42) is the collapse retardation factor, discussed in the penultimate paragraph 
of § 2.3. It follows from equation (|49p that, in the limit A ^ Aj,th, Tg, ~^ ^h,o- The growth time of this ambipolar- diffusion- 
induced fragmentation is shown in Figure 6c (the part of the curve labeled by "AD,fr"). Dimensionally, the growth time for this 
mode is VfffiT^^ ^ — r^^. ^/rnifi- It is the same as the nonZjnear solution found analytically by Mouschovias (1979; see also 1983, 
1987a, b, 1989) for the timescale of formation and evolution of protostellar cores (due to gravitationally-driven ambipolar 
diffusion). Numerical simulations (including the effects of grains, UV ionisation, and magnetic braking) of the formation of 
protostellar cores in magnetically supported molecular clouds have also found that the evolution occurs on this timescale 
(Fiedler & Mouschovias 1993; Ciolek & Mouschovias 1994, 1995; Basu & Mouschovias 1994, 1995a, b). The same timescale 
was found in the one-dimensional similarity solution of Scott (1984). It is clear from Figure 6c that, as predicted, f^^ would 
tend to 1/fni.o = 4.42 for A ^ A,j_th; see the inflection point in the curve (labeled "AD,fr"). However, f^^ falls below this 
would-be asymptotic value because, for A > Aj,mag, where Aj_mag is the magnetic Jeans wavelength (= 25.6; see eq. [54] below), 
gravitational forces on the neutrals, transmitted to the ions by neutral-ion collisions, exceed the restoring magnetic forces on 
the ions, and the ions and magnetic field are no longer able to remain stationary; the mode becomes a gravitational (Jeans) 
instability against the magnetic fleld, as originally found by Chandrasekhar & Fermi (1953); see Figure 6c, part of curve 
beyond the inflecton point, labeled "ff-|-". The approximations Vi^x — 0, — (which were used in deriving eq. [44 j ) are no 
longer valid, and equation (|49p no longer describes this mode. The proper expressions are derived below. 

The second mode resulting from the bifurcation of the modifled sound waves at As,n is described by the negative root of 
equation (|44|l . This is a neutral coUisional-decay mode (curve labeled "n,coH" in Fig. 66), with damping time 

2fni,0 



1+1 



Tni.O 



Aj,th/A)' 



1/2 







l--n% 


(A,T,th/A)'-l 



(50a) 



(50b) 



In equation (|50b|) we have again used the fact that, for the typical model cloud, f^; g <C 1. 

For A ^ Aj,th, "^d ■'ni.o. This limit is never attained by this mode, however, because the coUisional coupling between 
the neutrals and the ions becomes more effective with increasing A, and the ions (and, hence, the magnetic field lines) begin 
to move with the neutrals. As a result, the neutral coUisional-decay mode and the (ion) ambipolar-diffusion mode combine 
and merge to form wave modes, in a way similar to that for the neutral Alfven waves in § 3.1.2 (compare Figs. 46 and 66). 
This mode coupling occurs at A = 2.72 (= Ams.n; see below), which is only slightly greater than the neutral acoustic-wave 
cutoff Aa,n = 2.59 for the typical model cloud (see Fig. 66). The solution of the dispersion relation (eq. [IS]) for these modes 
in the limit |a;| <C l/fin,o is 



Cj — ±v. 



ms,n,0 1" 



16i^L,n,o 



1/2 



* -2 - r2 



where 



,1/2 



(51) 



(52) 



«ms,n,0 = (^A,n,0 + l) 

[or, in dimensional form, i;ms,n,o = {^A.n.o + C'a,o}^''^] is the magnetosonic speed in the neutrals. For the representative model 



* The nonlinear equivalent of this instant in the development of the ambipolar-diffusion-induced, gravitationally-driven fragmentation 
is the instant at which a fragment's mass-to- flux ratio reaches its critical value and dynamical contraction ensues with the magnetic field 
essentially frozen in the matter. 
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cloud used in this paper, iima.n.o 
< A < Aj,niag, where 



4.07 = l.OSiiA.n.o- Examination of equation (|51|) reveals that cDr > (i.e., > 0) if 



^«A,n,0^ni,0 



WA,n,0 



Aa,i 



''^ms,n,0 \'^ms,n.O^ 

is the neutral magnetosonic cutoff wavelength, and 

A.J. mag — STTl'ms.n.O 



(53) 



(54) 



(or, in dimensional form, A,j,mag = 27riims,n,oT"ft q) is the magnetic Jeans wavelength. For the typical model cloud, Ams.n = 2.72 
and Aj,mag = 25.6 (hence, Ams,n = 0.264 pc, and Aj,mag = 2.49 pc). [Note: in deriving equations ((53)) and ((54]), we have used 
the fact that, for the conditions of interest, fj^i g ^ 1 and l/i3ms.n,o ^ 1-] Inserting equations (|53[) and (|54p in equation (|5ip . 
the phase velocity for magnetosonic waves in the neutrals is found to be 

1/2 

, Ams,n A < Aj^mag (^^) 



1 



A 



A,, 



A„ 



A 



(see Fig. 6a, curve labeled "n,ms"). The waves are weakly damped by ambipolar diffusion (see Fig. 66, curve labeled "n,ms"); 
= 2AV47r^Oa (see eqs. [34b| and [5T]'). 

For A > Aj,mag, gravitational forces on the neutrals overwhelm the restoring magnetic forces, and the neutral magnetosonic 
modes become gravitationally suppressed (v^ = 0), just like the thermal Jeans modes become suppressed for A > Aj,th (see 
discussion in § 3.1.1). For A ~ Aj^mag, each (neutral) magnetosonic mode bifurcates into an ambipolar-diffusion mode (with 
damping timescale Tj — A^/47r^Oa for A 2> A.j.mag) and a conjugate Jeans mode (with damping timescale = 1 for 
A S> Aj,mag); see Figure 66, curves labeled "AD" and "ff— ", respectively. 

The conjugate Jeans (or classical cosmological) mode was discussed earlier in connection with equation psp . In this mode 
the neutrals, ions, and magnetic field lines are well coupled, but the gravitational forces prevent them from oscillating; motions 
are damped monotonically. The damping time, corresponding to the negative root of equation (|51|) . for A > A,T,mag is given by 

\ -1 



1 



J, mag 

A 



+ 



27r^Oa 
A2 



1/2 



27r^Pa 
A2 



(56) 



It is clear that, in the limit A — >■ oo, — >■ 1 (i.e., rd = % q); see Figure 66, curve labeled "ff— ". 

Finally, the gravitational instability mode (see eq. (49)) also changes behaviour at Aj,mag, as discussed above. For A > 
Aj,mag, the neutrals, plasma, and magnetic field lines are well coupled and behave like a single fiuid; self-gravity overwhelms 
restoring magnetic and thermal-pressure forces, and the mode behaves as a classical Jeans instability, in which density 
perturbations grow exponentially in time with a timescale 
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27r^Oa 



A2 



(57) 



which is equal to the damping time of the conjugate Jeans mode (see eq. [56]). From this equation we see that t^^ — >■ 1 (i.e., 
Tgr = ,-,) for A — >■ oo, in agreement with the long- wavelength behaviour of t^^ exhibited in Figure 6c. 



3.2.2 Longitudinal Modes for k _L Bq: Eigenvectors 
The main features of the eigenvectors are as follows. 

(a) The dominant component of the ion magnetosonic and the ion ambipolar-diffusion modes is (see Fig. 76, curves 
labeled "i,ms" and "AD"). It does not vanish at A = Aa.i; in fact, it hardly changes from 1, because the ion-AD mode 
maintains Vi^x large beyond Aa.i. Only when the ion-AD mode induces motions in the neutrals does Vi^^ begin to decrease, 
as the magnitude of Vn,x increases as A — >■ Ams.n (see Figs. 76 and 7a, curves labeled "i,ms" and "AD"). The ^-component of 
the magnetic field Bz in Figure 7c actually does not vanish. It is equal to Vi^x/vi^ms (see eq. (35]), but ■Ci,ms — 'i'A,i,o = 2.81 x 10'^. 

(b) The velocity Vi^x is also the dominant component of the eigenvector of the ion coUisional-decay mode at all wavelengths 
(see Fig. 76, curve labeled "i,coir'). 

(c) The eigenvector of the neutral acoustic mode has significant components iin,x and pn- At As.n, beyond which neutral sound 
waves do not exist and at which the acoustic mode bifurcates into the neutral collisional-decay mode and the pressure-diffusion 
mode (see Fig. 66), the collisional-decay mode is responsible for the increase in \vn,x\ as |pn| decreases to zero. Because motion 
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is induced in the ions as A approaches As.n (see Fig. 7b, curve labeled "n,coU"), \vn.x\ does not reach unity. 

(d) The most significant component of the PD mode is pn (see Fig. 7d); pn increases as A increases from As,n to Aj_th, at 
which wavelength pn reaches a maximum. Beyond Aj,th, ambipolar-diffusion-induced fragmentation sets in. At exactly Aj.th, 
Pn is the only nonvanishing component of the eigenvector of the AD,fr mode (see Figs. 7a - 7d). As A increases beyond Aj,th, 
the field lines begin to be compressed as the neutrals begin to couple to the ions, and Bz increases (see Fig. 7c) - at the 
expense of pn - while Vi^^ and Vn,x are negligible. As A approaches Aj^mag, the ambipolar-diffusion-fragmentation mode induces 
significant velocities Vn,x and Vi^x, while pn remains large (see Figs. 7a, 7b and 7d). As discussed in relation to Figure 6c, 
the AD,fr mode turns into the Jeans free-fall mode (ff-|-) beyond A ~ Aj,mag, as is clearly shown in Figures 7a and 7b (see 
curves labeled "AD,fr" and "ff+"). The bifurcation of the neutral magnetosonic mode into the AD and conjugate Jeans (or, 
cosmological, ff-) modes beyond A ~ Aj^niag, discussed in relation to Figure 66, is also seen clearly in Figures 7a and 7b. [Note: 
If we had plotted only the real part of the eigenvector, we would have found, for the neutral magnetosonic mode, that Re{pn} 
decreases discontinuously to zero, then increases smoothly, reaches a maximum, and then vanishes again as the magnetosonic 
waves are suppressed by gravity.] 

(e) Slightly beyond Ams,n, magnetosonic waves exist in the neutrals but, in addition to maintaining v-n,x large, they cause 
significant motion in the ions as well (see Figs. 7a and 76, curves labeled "n,ms"). The quantities Bz and pn are also 
nonnegligible. Note that at Aj,th, v-n,x and pn decrease while Vi^x increases. 

3.2.3 Transverse Modes for k ± Bo 

Examination of equation (|41cp reveals that there is one trivial mode with to = 0. The remaining equations governing the 
transverse modes with motions in the y-direction (eqs. [41a) and [41b] ) are identical with those in the ^-direction (eqs. [42a] 
and [42b]). Solving this simple system, we find that 



both of which are independent of wavelength. (Note: there are four modes in all, two corresponding to the first solution, and 
two corresponding to the second.) In the first mode the ions and neutrals move together with v^.y = Vi.y (or Vn,z ~ ^i.z); hence, 
there are no frictional forces between the two species, and f^j = oo. The second mode consists of oppositely fiowing streams 
of ions and neutrals. The momentum of each species decays by collisions with the other; hence, ion-neutral and neutral-ion 
collisions occur in parallel, and the net damping time for this mode is the harmonic mean of fin,o and fni,o (see eq. [58b] ). 
For the typical model cloud, the ion fluid has a much smaller inertia than the neutral fluid and, therefore, the streaming 
velocity of the ions is much greater than that of the neutrals. As a result, the motion of the ions is collisionally damped on 
the timescale ~ fin.o- 

3.3 Propagation at Angles 0° < 6 < 90° with respect to Bo 

Equations (]10a[) - ()10k[) reveal that, for 0° < < 90°, the equations for the modes with motions in the y-direction are 
uncoupled from those with motions in the (x, z)-plane. Motions in the x- and jz-direction, however, are coupled. The modes 
with Vn,y 7^ 7^ Vi^y (governed by eqs. [lOf] - [lOh] ) are purely transverse; those with nonvanishing velocities in the (x, 2:)-plane 
(governed by eqs. [10a] - [lOe] and [lOi] - [10k] ) are neither purely longitudinal nor purely transverse. 

3.3.1 Transverse Modes 

Figures 8a and 86 display the magnitude of the phase velocities and the damping timescales of the three transverse modes 
for propagation at an angle 9 = 45° with respect to Bo. Because they are incompressible, none of these modes can become 
gravitationally unstable. Eigenvectors are displayed in Figure 9. 

The dispersion relation for ion Alfven waves, obtained from equations (]10f[) - (]10h[) . has solutions given by equation l]24p . 
with WA.i.o cos 9 replacing UA.i.o in that expression. Thus, the two ion Alfven waves propagate with 



cD = 0, 



(58a) 



and 




(58b) 




21 1/2 



(59) 



(see Fig. 8a, curve labeled "i,A"). For A/(AA,iCos^) ^ 1, 5^ ~ iuA.i.ocosS = ±1.99 x 10^ for the typical model cloud. The 
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Figure 8. Eigenvalues of transverse modes as functions of wavelength, at an angle of propagation 8 = 45° with respect to Bq. All 
normalisations and labels are as in Fig. 2 and Table 1. There are three different modes, (a) Absolute value of phase velocity, \v^\. Phase 
speeds (large black circles) obtained from eq. I I59II are shown, in addition to those calculated from eq. I I61I I (open circles with interior 
black circles). (6) Damping timescales t^. Also depicted are values (downward- facing triangles with x's inside) obtained from eq. I I6UI1 . 
and values (diamonds with interior crosses) obtained from eq. I I62II . 
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Figure 9. Magnitudes of eigenvectors of transverse modes for 6 = 45° as functions of wavelength. All labels and normalisations are as 
in Fig. 3 and Tabic 1. As in Fig. 8, there are three modes, (a) Transverse (y-) component of the neutral velocity, lnn.j/l- {b) Transverse 
{y-) component of the ion velocity, \vi^y\. (c, d) Transverse (j/-) component of the magnetic field, \By\, shown as in Figs. 5c and bd. 
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Figure 10. Eigenvalues of transverse modes as functions of wavelength, at an angle of propagation = 10° with respect to Bq- All 
normalisations, labels, and symbols are as in Fig. 8 and Table 1. There are three different modes, (a) Absolute value of phase velocity, 
(b) Damping timescales t^. 



28 Mouschovias, Ciolek, and Morton 




Figure 11. Eigenvalues of transverse modes as functions of wavelength, at an angle of propagation = 80° with respect to Bq- All 
normalisations, labels, and symbols are as in Fig. 8 and Table 1. There are three different modes, (a) Absolute value of phase velocity, 
(b) Damping timescales t^. 
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± sign denotes degeneracy with respect to the direction of propagation. The damping timescale for these waves is 2Tin,o (see 
Fig. 86, curve labeled "i,A" ) . 

Collisions with the neutrals cut off the propagation of these waves for A > Aa.i cos 6 = 1.11 X 10"^ (see Fig 8a). For 
wavelengths greater than this value, each mode bifurcates into an ion coUisional-decay mode and an (ion) ambipolar-diffusion 
mode. The damping timescale of the former mode is given by equation (|27|) and is shown in Figure 8b (curve labeled "i,A"); 
it is identical with of the normal Alfven waves propagating along Bo, except for the fact that Aa.icosS replaces Aa.i in 
that expression (see discussion following eq. [24]). For very large A, = fin,o (see Fig. 8&, curve labeled "i,coU"). For the 
ambipolar-diffusion mode (Fig. 86, curve labeled "AD"), the damping timescale is the same as that given by equation (|28() . 
but again Aa.icosS replaces Aa,!; hence, for long wavelengths, it follows from equation H29ap that 



f, ~ . (60) 

47r2DaCos2 6l 

At small wavelengths the third mode is a low-frequency neutral coUisional-decay mode, with vi^y ~ (see Fig. 96, curve 
labeled "n,coU"). The purely imaginary frequency of this mode is given by equation (|3Up . iD = —i/T^ifi', hence the damping 
timescale is = fni,o {= 0.226 for the typical model cloud; see Fig. 86, curve labeled "n,coU"). At longer wavelengths, the 
motions of the neutrals and the ions become better coupled, and the neutral coUisional-decay mode merges with the ion 
ambipolar-diffusion mode (see Fig. 86), as discussed in §§ 3.1.2 and § 3.2.1. For longer wavelengths, Alfven waves in the 
neutrals can be sustained, as seen in Figure 8a. The frequencies of these modes are given by equations (|31ap and (|31b|) . but 
with tiA.i.o cos 9 and VA,n,o cos 9 replacing va,\,o and WA.n.o, respectively. The lower cutoff wavelength, below which these waves 
cannot propagate, is AA,nCos6' = 1.98. They have phase velocity 

1/2 

(61) 



= ±i;A,n,0 008 6^ 

and damping timescale 



Aa.ii cos t 



A 



47r2Da cos2 



(62) 



As A — >• CO, jC^I — ^ 'yA,n,o cosS — 2.79. 

The components of the eigenvectors of these modes displayed in Figures 9a - 9d are qualitatively similar to those shown 
in Figures 5a - 5d (for transverse modes propagating along iJo), but with (small) quantitative decrease in the values of the 
critical wavelengths. We therefore do not discuss them further for economy of space. 

In Figures 10a and 106 we display IS^I and f^, respectively, for the three transverse modes propagating at an angle 9 — 10° 
with respect to the magnetic field. The qualitative solutions of the dispersion relations for these modes are the same as for 
those described above for propagation at an angle 9 = 45°. Quantitatively, the differences in Figures 10a - 106 and 8a - 86 
stem from the different values of the numerical factors cos 9 and sm9. Similarly, the magnitude of the phase velocities and the 
damping timescales for the transverse modes propagating at an angle 9 = 80° with respect to Bq are shown in Figures 11a 
and 116. They again differ from those for 9 — 45° because of the factors cos 9 and sin 6. Although qualitatively the modes at 
9 = 10° and 9 — 80° are the same, comparison of Figures 10 and 11 reveals that the quantitative differences in the magnitude 
of the phase velocities and the cutoff wavelengths are substantial. Thus the critical wavelengths that determine which modes 
can or cannot propagate in molecular clouds are expected to be very different depending on their direction of propagation with 
respect to the mean magnetic field Bo- 



3.3.2 Modes with Motions in the (x, z)-Plane 

Figures 12a, 126, and 12c exhibit, respectively, the magnitude of the phase velocity v^, the damping timescale f^, and the 
growth time f^, as functions of A for the specific case of propagation at ^ = 45° with respect to Bq for modes with motions 
in the {x, z)-plane. Eigenvectors are displayed in Figures 13a - 13g. There are 7 modes in all displayed in these figures. 

As in the preceding sections, one of the ion modes is a coUisional-decay mode (see Figs. 12 and 13, curves labeled "i,coU"), 
with the ions streaming through a fixed background of neutral particles {vn,x = "Cn.z = 0); for this case Bz ~ 0, and the ions 
move with Vi^x = Vi,z. The frequency is given by equation (|16p : hence, the damping timescale is — fin.o {~ 4.45 x 10^^ for 
the typical model cloud). This is the horizontal line labeled by "i,coU" in Figure 126. 

There also exist two high-frequency ion wave modes. These waves are ion fast modes. For these modes, v\^x/vi,z ~ 
— sm9/ cos9 — —1 ai 9 = 45° (see Figs. 13c and 13d, curves labeled "i,fast"; in these Figures the "i,fast" curves coincide with 
the "i,coU" curves). Since the ion Alfven speed and the magnetosonic speed (in this typical case) are essentially the same, the 
dispersion relation describing them is the same as equation ((24]); it is again the case that the waves are cut off at A > Aa.i 
{— 1.57 X 10~^ ; see Fig. 12a). For wavelengths > Xaa, each ion mode bifurcates into an ion coUisional-decay mode and an ion 
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Figure 12. Eigenvalues of longitudinal modes as functions of wavelength, at an angle of propagation 8 = 45° with respect to Bq- All 
normalisations and labels are as in Fig. 2 and Table 1. There are seven different modes, (a) Absolute value of phase velocity, \v^\. Boxes 
with interior six-pointed stars denote phase velocity values derived from eq. I I66II . Phase speeds calculated by using eq. I I70II are shown 
as left-facing triangles, and phase speeds obtained from eq. {TTJ are shown as boxes with interior crosses. (6) Damping timcscalcs t^. 
Also displayed are values obtained from eq. Il29all (downward-facing triangles) and from eq. II71I I (left-facing triangles with interior black 
circles). Values obtained from eq. I|34a|l are also displayed as crosses, (c) Growth timescale t^^. 
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ambipolar-diffusion mode (see Fig. 126). The decay timescales for these two modes are the same as those previously discussed 
for the cases 6* = 0° and 90° (see § 3.1.2 and § 3.2.1). 

The fourth mode is the ion dissociative-decay mode, discussed previously in §§ 3.1.1 and 3.2.1. In this mode, nonpropa- 
gating density perturbations in the ions decay by dissociative recombinations of molecular ions and electrons. The decay time 
is = (2pi,oam,dr)~^ {~ 1.28 X 10~* for the typical model cloud); this is the horizontal line labeled by "i,rec" in Figure 12&. 

At small wavelengths there are two low-frequency acoustic-wave modes in the neutrals and one neutral coUisional-decay 
mode. The dispersion relation for the acoustic modes, which are predominantly polarised in the i-direction, is most easily 
found by first finding m^x in terms of Vn,x and uj. In the limit \uj\ <C l/rin.o and ujBz ^ 0, one finds from equations (|10i() . P0j[ ), 
(llOkD. and (llOdI) that 



; COS 



OJ 4- {i/fni,o) 



uj + (icos^ 6'/fni,o) 



(63) 



For > l/Tni,o, the term in brackets in equation (|63p is essentially unity, and Vi^x — Vn,x cos^9. The factor cos^S is a measure 
of the opposition presented to the neutrals by the ions, which are attached to magnetic field lines. If = 0° the ions are 
effectively inertialess, and the neutrals sweep them up easily, so that Vi^x = Vn,x- However, if 6 — 90°, Vi^x — 0, i.e., the ions 
are held in place by the magnetic field as the neutrals move through them. In this case, the ions present the stiffest opposition 
to the neutral motion. For > l/rni,o, we insert equation (|63[) in equation (|10cp to find the solution 



±k i 1 



1 



1 + 



sin ''6' 

2rni.O 



1/2 



. sin^ e 
■ 

2fni,0 



[Note that for 9 — 0° or 90°, equation (|64p reduces to equations (jlSp or (I44p . respectively.] For 
where 



sin^e + (2fni,o)'"''^' 



1 + (2fni,o)2 

these modes are sound waves, with 



\ As,n 



1/2 



(64) 



(65a) 



(65b) 



(66) 



(see Fig. 12a) and = 2fni,o/sin^S — 0.904 for 6 — 45° (see horizontal line labeled "acoustic" in Fig. 126). AB,ni6) is the 
acoustic-wave angular (or, stiffness) parameter: for 9 — 90° equations (|65ap and (|66p yield our earlier result, that sound waves 
propagate only for A < As.n, while for 9 = 0° we recover our other earlier result that Aj^th is the upper cutoff wavelength for 
these waves. 

The sound-wave upper cutoff wavelength (|65ap is relevant only if the waves do not transition into neutral slow modes 
at larger wavelengths (see discussion below). It turns out that this depends on the angle 9 of propagation: the upper cutoff 
equation is applicable only to sound waves propagating at 9 > Smax (for the typical model cloud, Smax ~ 62.9° — see below). 
For sound waves propagating at such angles, there is a mode bifurcation at the cutoff H65ap . The dispersion relation H64p 
reveals that at larger A one of the modes becomes a neutral collisional- decay mode, with damping timescale — > fni.o/sin^ 6, 
and the other mode becomes a pressure-driven diffusion mode with damping timescale 



(AAj.th)'sin2 6l 



for A < Aj,t 



(67) 



Op[l-(A/Aj,th)2] 

This diffusion timescale is the same as that of equation (|47ap . except that it is multiplied by sin^ 9, which reflects the reduced 
effectiveness of neutral-ion collisions in slowing down the neutrals at angles 9 < 90° with respect to the magnetic field. 
At A = Aj.th the timescale (|67p becomes infinite. For A > Aj_th there is a neutral ambipolar-diffusion-mdueed gravitational 
fragmentation mode at the angle 9 (> Smax), having a growth time 



sin^ 9/fnifi 



vsfi sin 9 



(A.7,th/A)' 1- (A,T,th/A)' 



(68) 



The fragmentation time (I68p is equal to the growth time (|49ll multiplied by sin 9, again reflecting the reduced collisional 
resistance on the neutrals by magnetically-coupled ions when the angle of propagation is less than 90° with respect to the 
field direction. For A 2> Aj,th, equation (|68p yields f^^ — > sin'^S/fni,o ~ vsfisiv? 9. However, similar to what occurs in the 
case for 9 = 90° (see § 3.2.1 and Fig. 5c), this limit will not be attained because the approximation of stationary field lines 
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Figure 13. Magnitudes of eigenvectors of modes with motions in the (x,z)-plane for = 45° as functions of wavelength. All labels and 
normalisations are the same as in Fig. 3 and Table 1. As in Fig. 12, there are seven modes, (a) Longitudinal (a;-)component of the neutral 
velocity, liin.xl; (b) l^n.^l- (c) Longitudinal (a;-)component of the ion velocity, \vi^x\', (d-) 
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Figure 13. Cont. (e) Neutral density |pn|- (/) Ion density |pi|. (g) Magnetic field component \Bz\. 

breaks down when A becomes > Aj,mag- When this happens, magnetic forces are overwhelmed by self-gravitational forces and 
— >■ 1 at these larger wavelengths. 
From equation (|64p one finds that collisional damping of the motion of the neutrals by ions causes Iw^l to become less 
than unity for A near the value given by the right-hand side of equation (I65a|) . However, for < ^ < 90°, the sound waves are 
not always cut off at this wavelength. This is due to the fact that the bracketed term in the expression for the a;-component 
of the ion velocity (eq. [63]) is no longer essentially unity (because jcD] ~ l/Tni,o) at these wavelengths; thus, the dispersion 
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relation (|64|) is no longer valid. Physically, this is a result of the fact that the ions move readily with the neutrals in the 
a;-direction at these frequencies; this means that the waves suffer less damping, because the frictional force on the neutrals 
is reduced when the ions and the neutrals move together. (We note that, for 9 < 90°, equation [63] yields Vi^x — Vn.x in the 
limit \uj\ < l/fni,o.) 

For A sufficiently large, such that \Cj\ < cos^ 0/Tni,o, the frequency of the neutral waves is 



u ~ ±fc cos 6 



1-1 

fc2 
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T-aifik sin 6 tan 6 
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Hence, for these conditions, the phase velocity is 
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(71) 



These waves are neutral slow modes, modified by gravity, and cannot propagate at wavelengths A > Aj^th {= Stt for the typical 
model cloud); at ^ = 0° the slow mode dispersion relation (|69[) is identical to the relation for the Jeans mode (eq. [H]). From 
the low-frequency condition used to derive equation (I69p it is found that slow modes exist only for wavelengths 
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The quantity S-niO) is the angular slow- mode factor. The relation (|73p is derived from the requirement that the slow modes 
arise from the acoustic wave modes. For this mode conversion to occur, the minimum slow mode wavelength (I72a|) must be 
less than or equal to the acoustic- wave upper cutoff wavelength (|65ap : thus, the inequality H73[) follows. This is equivalent to 



having 



+ (2fni,0)2 



where ^max is defined by the condition 
1 

4 ' 



(74) 



\i 9 < ^max, slow modes emerge from the acoustic modes (without a bifurcation) and propagate for As,n5n(^)^A < Aj.th. 
Otherwise, when 6 > Smax, the acoustic waves have an upper cutoff and mode bifurcation occurs at the wavelength ((65a]); in 
that case, there are no slow modes. 

For the typical model cloud, with fni,o = 0.226, Smax = 62.9°. At ^ = 45° the slow mode minimum wavelength 
As,n5n(45°) = 2.01, and the transition from sound waves to slow modes can be seen in Figures 12a and 126 (curves la- 



beled "n,slow") to occur at this wavelength; Figures 13a and 13c show that, for these modes, 



at wavelengths 



greater than or equal to the transition wavelength. 

Note that, as A -> Aj^th from below, equation shows that the slow mode damping timescale f ^ — >■ co (see curve 
labeled by "n,slow" in Fig. 126). For A > Aj_th there are again the two conjugate modes: the gravitational instability (or 
Jeans) mode (see Fig. 12c, curve labeled "ff-|-") and the "cosmological" mode (see curve labeled "ff— " in Fig. 126; although 
this curve "crosses" the curve labeled "n,fast" in Fig. 126, the two modes do not actually interact); as A — >■ co, both the 
growth timescale of the unstable mode and the damping timescale of the cosmological mode go to unity (i.e., in dimensional 
form, Tgr = Td = Tjf o), as seen in Figures 126 and 12c, respectively. 

The other mode affecting the neutrals at short wavelengths is a neutral collistonal- decay mode ("n,coU"). The velocity 
of the neutrals at small A is predominantly in the z-direction for this mode (see Figs. 13a and 136); the frequency is purely 
imaginary and given by 

ijj = —I— . (75) 



For clouds with perfect neutral-ion coUisional coupling, i.e., f^\fi = 0, eq. II74I I yields Sn 



: 65.5° 
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Figure 14. Eigenvalues of longitudinal modes as functions of wavelength, at an angle of propagation 6 = 10° with respect to Bq- All 
normalisations, labels, and symbols are as in Fig. 12 and Table 1. There are seven diflferent modes, (a) Absolute value of phase velocity, 
\v^\. (b) Damping timescales Tj. (c) Growth timescale Tg,. 
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Figure 15. Eigenvalues of longitudinal modes as functions of wavelength, at an angle of propagation 8 = 80° with respect to Bq. All 
normalisations, labels, and symbols are as in Fig. 12 and Table 1. There are seven different modes, (a) Absolute value of phase velocity, 
\v^\. (b) Damping timescales . Five-pointed stars with interior black circles represent values calculated by using eq. II67I I . ( c) Growth 
timescale t . Values obtained from eq. II68I I are depicted as six-pointed stars with interior black circles. 
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Thus, = fni.o/cos^ e = 0.452 (see Fig. 126, curve labeled "n,coU"). For A greater than the value given by the right-hand 
side of equation (|65ap . it is again the case that the motion of the ions and magnetic field becomes better coupled to that 
of the neutrals. In this wavelength regime \vn,x\ — l^n.z] — — |Ci,z| (see Figs. 13a - 13rf). This mode merges with the 

ion ambipolar-diffusion mode at A = Ams,n (see Fig. 12b), and, for A > Ams.n, fast modes are excited in the neutrals (curves 
labeled "n,fast" in Figs. 12a, 126, and 13a - 13d), degenerate with respect to the direction of propagation. In these modes, 
the polarisation is given by Vn,x/vn,z = — sinO/ cosO (= —1 for 9 = 45°). Hence, Vn ■ Bo — 0, i.e., the fast modes are polarised 
perpendicular to the magnetic field. The dispersion relation for these modes is essentially the same as equation (|51|) because 
of the fact that VA,n o 3> 1. They decay, as they propagate, on the ambipolar-diffusion timescale (see Fig. 126). For A > Aj,mag, 
these waves tend to get suppressed by gravity. 

For A > Aj,mag thermal and magnetic restoring forces in the a;-direction are overwhelmed by gravitational forces, making 
Vn,x — (see Fig. 13a). Hence the modes become essentially incompressible. Waves are still able to propagate for longer 
wavelengths, however, because of the transverse restoring magnetic tension force (i.e., ^ 0; see Fig. 13g). Solving equations 
(|10a|) - (|10d[) and (|10i|) - (|10k|) with the conditions Vn.x = — pn in the limit \uj\ <^ l/rin.o, we find 
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Hence, these modes are modified Alfven waves in the neutrals, with 
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2.79, in agreement with the long-wavelength behaviour of 
n,A"). This is yet another example of a transition between 



and = 2A^/47r'^I?a. In the limit A — >■ oo, \v^\ ~ ?}A,n,ocosS 
these modes shown in Figure 12a (curves labeled "n,fast" and 
wave modes without bifurcation. 

Figures 14a, 146, and 14c show IC^I, f^, and fg, as functions of A for the seven different modes with motions in the (x, z)- 
plane propagating at an angle 6 = 10° with respect to Bq. Comparison with Figures 12a - 12c reveals that the qualitative 
behaviour of the various modes as functions of wavelength is the same as in the case of propagation at S = 45° . The quantitative 
differences stem from the numerical factors cos 6 and sin 6, which become substantial for 6 approaching 0° or 90°. As Figure 
14a shows clearly, wave modes in the neutrals exist at all wavelengths and their decay times are very long (see Fig. 146, 
curves labeled "acoustic", "n,slow", and "n,fast"). Figures 15a, - 15c show the same quantities as Figures 14a - 14c but for 
propagation at = 80° with respect to the unperturbed magnetic field Bo. There are no slow modes in Figure 15 (unlike 
the cases in Figs. 12 and 14) because 6 > Smax for that angle of propagation in the typical model. Instead, the sound waves 
are cut off at the maximum wavelength A = As,n/.4s,n(80°) = 2.66, where there is a bifurcation. At wavelengths greater than 
this maximum, the modes are a pressure-driven diffusion mode ("PD" curve) and a neutral coUisional-decay mode ("n,coU"). 
There is also an ambipolar diffusion-induced fragmentation mode seen in Figure 15c ("AD,fr" curve), which approaches the 
predicted limiting value (see eq. [68]) of sin^ 80°/fni,o ~ 4.29 at A just below Aj^mag (= 25.6 in the typical model cloud). 



4 SUMMARY AND DISCUSSION 

We have obtained and analyzed the dispersion relations for MHD wave modes and instabilities for different directions of 
propagation with respect to the zeroth-order magnetic field Bo in a two-fluid weakly ionised system, and we have applied 
the results to a typical interstellar molecular cloud. The system of equations has four dimensionless free parameters, fni.o, 
"^in.o, CA,i,o, and fim.dr. They represent, respectively, the neutral-ion (momentum-exchange) collision time and the ion-neutral 
collision time in units of the free-fall time of the zeroth-order state, the Alfven speed in the ions in units of the adiabatic speed 
of sound in the neutrals, and the dissociative recombination coefficient (see eq. [lld| ). (Because of ionisation equilibrium in 
the zeroth-order state, the dimensionless cosmic-ray ionisation rate C,^-g_ is expressible in terms of (3m,dr, ■fln.o and Tni.o.) 

There are two distinct kinds of ambipolar diffusion, whose combined effect is unavoidable in typical molecular clouds and 
has crucial consequences on their evolution: 

(a) In the presence of hydromagnetic (HM) waves or turbulence, the tension of field lines (or the outward pressure due to com- 
pressed field lines) drives the motion of charged particles relative to the neutrals, with the tendency /consequence to straighten 
out the bent or tangled magnetic field lines (or to move compressed field lines apart, toward a more uniform configuration). 
The timescale of this process is proportional to the square of the wavelength of the HM waves (or the characteristic length of 
the field-line tangling, or the magnitude of the field gradient) - see eqs. (|34ap and (|80p . For lengthscales typical of molecular 
cloud cores (< 0.1 pc), it is much smaller than the free-fall time. This is the magnettcally-driven ambipolar diffusion. It is 
this kind of ambipolar diffusion which is responsible for the observed large-scale ordering of polarisation vectors, indicating 
large-scale ordering of the magnetic field lines in molecular clouds. 
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(b) Gravitationally- driven ambtpolar diffusion sets in with a short enough timescale, but longer than the free-fall time, in 
the deep interiors of self-gravitating clouds, where the degree of ionisation is Xi < 10"'^. Its onset can be spontaneous or 
initiated as a result of magnetically-driven ambipolar diffusion depriving a sclf-gravitating cloud of any support that most 
short- wavelength HM waves (or turbulence) may have provided against the cloud's self-gravity (Mouschovias 1987a). It (i) 
allows the clouds to fragment as neutral particles contract through almost stationary field lines (and the attached charged 
particles); and, consequently, (ii) increases the mass-to-flux ratio of the forming fragments (or cores). Once the critical mass- 
to-flux ratio (Mouschovias & Spitzer 1976) of fragments is exceeded, dynamical contraction ensues, while the cloud envelopes 
remain magnetically supported, as found by numerical simulations starting with Fiedler & Mouschovias (1993) and confirmed 
by numerous observations. 

Hydromagnetic (HM) waves with phase velocity 

can propagate in all directions with respect to Bo, provided that A > Aa.h, where 

Bo \ /2 X 10^ cm-^\^''V2 X 10"' 



The long-wavelength waves are long-lived; the decay time is the magnetically-driven ambipolar-diflusion timescale 

which is to be distinguished from the growth time of gravitationally-driven ambipolar-diffusion, relevant to fragmentation of 
molecular clouds into self-gravitating cores; namely, 

rAD,f. = ^^l.lxlO« (-^) yr. (81) 

The (one-dimensional) free-fall timescale at the density nn,o = 2 x 10'^ cm^'^ as A — >• cxa is % ^ = (47rGpn,o)"^^^ = 3.9 X 10^ 
yr. The nonlinear counterparts of these modes have been shown to explain quantitatively the observed highly supersonic but 
subAlfvenic linewidths in molecular clouds, their cores, and even in OH and H2O masers in which the strength of the magnetic 
field has been measured (Mouschovias & Psaltis 1995; Mouschovias et al. 2006). 

Most HM waves with A < AA,n cannot propagate in the neutrals because they are damped rapidly by ambipolar diffusion. 
This means that there cannot be any contribution from this wavelength regime to the spectrum of HM "turbulence" in 
molecular clouds (Mouschovias & Psaltis 2011, in preparation), which may provide clouds with a source of nonthermal 
pressure. This led Mouschovias (1987a) to argue that the decay of HM waves by ambipolar diffusion on wavelength scales 
< 0.1 — 0.3 pc can initiate the formation of protostellar cores in otherwise magnetically supported clouds (see also Mouschovias 
1991a). Damping of short-wavelength HM waves by ambipolar diffusion has also been proposed (Mouschovias 1987a, § 2.2.5) as 
the cause of the observed narrowing and thermalization of linewidths with increasing column density, as observed, for example, 
by Baudry et al. (1981) in the cloud TMC 2. Depending on the angle of propagation with respect to the unperturbed magnetic 
field Bq, however, we find that certain long-lived modes in the neutrals exist at all wavelengths, while ion modes usually damp 
very rapidly even at short wavelengths. This may explain the observations by Li & Houde (2008) in the M17 molecular cloud, 
which show neutral motions at small lengthscales but much smaller ion motions on the same scales. 

Gravitational instability is found to set in at A = Aj^th (the thermal Jeans wavelength) for all 9, in agreement with 
Chandrasekhar & Fermi (1953). For 9 ~ 90°, though, the timescale for the instability is ~ t'ff.oTff^o = /''"">. O' which is 
the gravitationally-driven ambipolar-diffusion timescale. This is in agreement with the result for the ambipolar-diffusion- 
initiated formation and contraction of protostellar fragments (or cores) obtained from nonlinear calculations analytically by 
Mouschovias (1979, 1989) and Mouschovias & Paleologou (1986), and numerically by Fiedler & Mouschovias (1993) for a 
two-fluid system (neutral molecules and plasma), Ciolek & Mouschovias (1994) for a four-fluid system (i.e., neutral molecules, 
plasma, negatively-charged, and neutral grains) , and Basu & Mouschovias (1994, 1995a, b) for a two-fluid system inchiding 
rotation and magnetic braking. It is also in agreement with the results of simulations by Tassis & Mouschovias (2007a, b, 
c) and Kunz & Mouschovias (2009, 2010) for a six-fluid system (neutral molecules, atomic and molecular ions, electrons, 
negatively-charged, positively-charged, and neutral grains). 

Another recent success of the linear theory (applied to model molecular clouds flattened along the magnetic field) in 
predicting inherently nonlinear phenomena was demonstrated by Kunz & Mouschovias (2010). They obtained analytically 
the Core Mass Function (CMF) resulting from gravitationally-driven ambipolar-diffusion-induced fragmentation of molecular 
clouds. The predicted CMF is in excellent agreement with observations of more than 300 cores in Orion (Nutter & Ward- 
Thompson 2007), not only at the high-mass end (for which many alternative models can obtain such agreement), but also at 
the turnover mass and low-mass end. 
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Table 2. Critical Wavelengths 



Dimensionless Expression 



Model Valuet 



Dimensional Expression 



ion Alfvcn-wavc upper cutoff 

ion magnctosonic-wavc upper cutoff 

acoustic-wave upper cutoff 

neutral Alfven-wave lower cutoff 

neutral magnetosonic-wave lower cutoff 

thermal Jeans 

magnetic Jeans 



AA,i = 47r'DA,i,0'nn,0 
AA,i 

As,n = 47rfni,o/[l + {2fni,o)^ 
AA,n = T^VA,-n,OT-n\fi 
Ams.n = (i'A,n,o/i'ms,n,o)AA. 
Aj,th = 27r 
Aj^mag — 27ri^ms,n,0 




0.0157 
0.0157 
2.59 
2.80 
2.72 
27r 
25.6 




^To convert to parsecs, multiply by the unit of length for the typical model cloud, Ca,OTfi,o = 9-72 X 10 ^ pc. 

Table 2 lists all the critical wavelengths present in a two-fluid system, such as a typical molecular cloud. The name of each 
critical wavelength, its defining dimensionless expression and typical value, and the corresponding dimensional expression are 
listed in columns 1-4, respectively. 

Table 3 summarizes conveniently all the modes that can exist in a two-fluid system (typical molecular cloud) for propa- 
gation parallel, perpendicular, and at an arbitrary angle with respect to the unperturbed magnetic field. The name of each 
mode is shown in the first column, and the wavelength range in which the mode exists is shown in the second column. 

Although in a single-fluid system linear modes are independent of one another, in a multifluid system (such as a molecular 
cloud) this is not the case. A mode in one fluid can bifurcate due to interaction with the other fluid and give rise to two 
distinct daughter modes. Such interaction is also responsible for the reverse phenomenon of mode merging. Moreover, in a 
multifluid system, a wave mode in one fluid/species (e.g. a sound wave in the neutrals) can transition into another wave 
mode (e.g., a slow MHD wave in the neutrals) without bifurcation. In other words, linear waves in a multifluid system exhibit 
behaviour that only nonlinear waves exhibit in a single-fluid system. 

In a following paper we present a study of the free parameters, spanning the range of observationally allowed values, to 
examine their effect on HM waves and instabilities in molecular clouds. 

This work was carried out in 1987 - 1989, when the bulk of this paper was also written. The paper was updated in 
1994 - 1996, but never submitted for publication due to spatial dispersal of its three authors and their involvement in code 
development for nonlinear calculations. Its Introduction and some other parts have been updated to refer to relevant work 
published since then. The work was supported in part by the National Science Foundation under grant AST-93-20250, and 
also by the New York Center for Astrobiology under grant # NNA09DA80A. TChM acknowledges the hospitality of the 
Alexander von Humboldt Foundation and of the Max-Planck-Institut fiir Extraterrestrische Physik during the writing of part 
of this paper in 1994. GC is supported by the New York Center for Astrobiology (a member of the NASA Astrobiology 
Institute) under grant # NNA09DA80A. 
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